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Abstract 



This thesis examines three problems of the quark model, split between two general 
areas. 

The first area involves models of meson decay. The calculations of decay widths in 
two such models, the ^Pq model and the flux-tube breaking model, were automated 
so that they could be easily done for any decay. The models were then used to 
investigate two particular problems. The first is the nature of the /4(2220) - although 
tentatively identified as the ss meson by the Particle Data Group, its identity has 
been uncertain since its discovery in 1983. An exhaustive calculation of the strong 
decay modes of the and ^F^ ss mesons was performed in order to examine the 
possibility that the state is a meson. It was found that the /4(2220) cannot be the 
■^^2 SS meson, and is unlikely to be the "^^4, although given the uncertainties of the 
models the latter possibility cannot be ruled out. Instead the following explanation is 
proposed: that the broad state seen in hadron beam experiments is the ss meson, 
and the narrow state seen in J/ radiative decay is a glueball. Further experimental 
data is needed to finally identify the /4(2220). 

The second problem investigated using meson decay models is the determination 
of the mixing angle between the i^ri(1270) and i^ri(1400) mesons. This was done by 
comparing predictions of the meson decay models for five partial decay widths and two 
ratios of D to S amplitudes to experimental data. A mixing angle of approximately 44° 
was extracted, and the implications to the quark model Hamiltonian are discussed. 

The third problem falls in the second general area: final state interactions (FSl's). 
The effect of the strong interactions between the outgoing pions was examined for 
the reaction 77 — > tttt near threshold. For charged pions, the experimental data 
agrees with the results calculated both with and without the effects of FSI's - better 

iii 



data is need to distinguish between the two. However, agreement was not obtained 
for neutral pions. It is beheved that the discrepancy may be due to the effects of 
resonance production. If this proves to be the case, then the techniques used in this 
work could be extended to try to understand a long-standing puzzle: the nature of 
the unexplained structures seen in the cross sections of 77 — > two vector mesons. 
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Chapter 1 
Introduction 



Particle physics is the study of the particles that make up the universe, and the 
interactions that take place between them. It is also referred to as elementary particle 
physics, by which we mean that we would ideally like to understand the universe 
in terms of its "elementary" constituents (should such things exist), that are not 
composed of other particles. 

Why study particle physics? The questions "What is the universe made of?" 
and "How does it work?" are of fundamental interest for their own sake. In addition, 
understanding the universe as it is now might help us answer the question "Where 
did it come from?" 

Aside from these weighty questions, what about immediate material benefits? 
Like much of fundamental research, it is impossible to know now what benefits might 
be realized from particle physics. However, it is worth noting that in 1897 the quest 
to understand the universe led to the discovery of the electron by J.J. Thompson. 
It is safe to say that this early particle physicist could not have predicted the huge 
effect his work would have on humanity. 
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1.1 The Standard Model of Particle Physics 

We assume that the reader is famihar with modern physics. Famiharity with the field 
of particle physics would also be useful, but hopefully not necessary; in addition to 
the introduction to the field that follows, Appendix ^ contains descriptions of some 
of the terms and tools of particle physics that may be useful to specialists in other 
fields. 



1.1.1 The Particles that Constitute Matter 

At present, the elementary particles that make up matter are thought to be the 
quarks (six of them, plus their antiparticles) and the leptons (again, six of them, plus 
antiparticles) . We think of them as elementary because, so far, we do not have any 
indication that they are not. The quarks are the up (u), charm (c) and top (t), all 
with electric charge +| (on a scale where the electron has charge —1), and the down 
(d), strange (s) and bottom (6), all with charge — |. The leptons are the electron (e), 
mu (/i) and tau (r), all with charge —1, and the three corresponding neutrinos, u^, 
and all with charge 0. We say corresponding because all these particles are 
usually thought of as being arranged in doublets. 



Quarks : 



charge + 
charge — 



2 . 

3 ■ 

1 . 

3 ■ 




Leptons 



charge : 



/ 



T 



charge — 1 : \ ^ j \ 

As one moves to the right in the first, second and fourth rows^, the particles become 
more massive. The particles in the third row, the neutrinos, have masses smaller 



^Each column is called a generation; there appear to be only three generations. 
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Quark 


Quark Mass 


Lepton 


Lepton Mass 




Current 


Constituent 






u 


2 to 8 MeV 


~ 300 MeV 




< 5.1 eV (95% CL) 


c 


1.0 to 1.6 GeV 


~ 1.5 GeV 




< 0.16 MeV (90% CL) 


t 


180 ± 12 GeV 


~ 180 GeV 




< 31 MeV (95% CL) 


d 


5 to 15 MeV 


~ 300 MeV 


e 


0.51099906(15) MeV 


s 


100 to 300 MeV 


~ 500 MeV 




105.658389(34) MeV 


b 


4.1 to 4.5 GeV 


~ 5 GeV 


T 


mr.ltlt MeV 



Table 1.1: Masses of the quarks and leptons. The numbers in parentheses denote the 
1 standard deviation uncertainty in the end digits of the value. The ± errors are also 
1 standard deviation uncertainties. CL means confidence limit (e.g. we are 95% sure 
that rui,^ < 5.1 eV). 

than the experiments can currently measure, and may be massless. The present 
experimental values of the masses^ are listed in Table |1.1| ||T|. Two different masses 



are listed for the quarks. The current mass is the mass that appears in the Lagrangian 
describing the strong interaction (see below). The constituent mass is the effective 
mass the quark has when it is bound inside a hadron (the hadrons are the strongly 
interacting particles - see below) - the numbers given are approximate because they 
depend on the hadron model used. 

Associated with each type of quark and lepton is its antiparticle, which has 
the same mass, but opposite quantum numbers (such as charge). Antiparticles are 
denoted by an overbar (e.g. u). Each type of quark or lepton is referred to as a 
flavour; a particle and its antiparticle have opposite flavour. Each quark and lepton 
has an intrinsic angular momentum, called spin, of |, making them fermions (they 
obey Fermi-Dirac statistics, which hold for particles of half-integral spin). 

In addition to their electric charge, each quark has an additional "charge" referred 



^Note that the masses in Table 1.1 are given in units of energy (eV), when the units of mass would 
normally be eV/c^. We are following the particle physics convention of setting c = h = 1, which 
expresses mass and momentum in units of energy, and length and time in units of energy"^. When 
converting from these 'natural units' to more standard units at the end of a calculation, factors of 
Jt and c are inserted as needed to give the required units. 
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to as colour (but having absolutely nothing to do with the colours of the everyday 
world). There are three possible values of a colour charge, plus the three anti-colours 
of the antiquarks. It appears to be a property of nature that coloured objects cannot 
exist freely by themselves, so quarks are confined inside hadrons in configurations 
that produce an object with no net colour. 

1.1.2 The Interactions 

There are four different types of interactions that take place between the particles; to 
each of them there corresponds one or more gauge bosons (they obey Bose-Einstein 
statistics, which hold for particles of integral spin), a type of particle that carries the 
effects of the interaction when it is exchanged between the two particles involved. 

The electromagnetic interaction occurs between objects with electric charge, in- 
cluding neutral particles composed of charged constituents (e.g. the neutron interacts 
via its magnetic moment). It is carried by the photon 7, which is massless and has 
a spin of 1 (making it a vector particle). As a quantum of electromagnetic radia- 
tion, it will be most familiar to the reader as visible light. The theory governing the 
electromagnetic interaction is called Quantum Electrodynamics (QED). 

The weak interaction occurs between leptons and quarks, and involves the in- 
teracting particles being transformed into their partners in the doublets of particles 
shown in ( |1.1| ).'^ It is carried by the and bosons (spin 1) with the charges 
shown and masses of 91.187±0.007 GeV and 80.33±0.15 GeV, respectively. Because 
the weak interaction deals with the flavour of the involved particles, the theory that 
describes it is sometimes called Quantum Flavourdynamics. The theories of the elec- 
tromagnetic and weak interactions have been combined into a single theory, called 

•^For the quarks at least this is a simphfication: while the quarks prefer to be transformed into 
their partners, they may also be transformed into one of the other two quarks in the same row as 
their partners. It is not clear if the leptons also have this type of cross-generational mixing. 
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the Electroweak Theory. The weak interaction is the cause of certain decays, such as 
the decay of the neutron n — > p e~ Pe- 

The strong interaction occurs between coloured objects, including colourless par- 
ticles composed of coloured constituents. It is carried by eight gluons g, which have 
spin 1, no mass and no charge. However, they themselves have colour, which compli- 
cates the theory of the strong interaction. Quantum Chromodynamics (QCD). The 
strong interaction binds the quarks together to form hadrons, including the nucleons. 
The residual strong forces between the colourless nucleons bind the nucleus together. 

The particles and interactions discussed above, and a spontaneous symmetry 
breaking mechanism (called the Higgs mechanism) that gives the fermions their 
masses make up what is called the Standard Model of Particle Physics. The Higgs 
mechanism gives rise to the Higgs boson, a massive (> 58.4 GeV, 95% CL), neutral 
scalar (spin 0) particle, which has not yet been seen. 

The remaining interaction, gravity, is well described by a classical theory, General 
Relativity, but so far no way has been found to incorporate it into a quantum field 
theory (such as the electroweak theory or QCD). That has not prevented physicists 
from naming its hypothetical carrier the graviton - it would be massless, neutral, and 
have spin 2. 

One drawback of the standard model is that it has 19 (23 if the neutrinos have 
mass) parameters that must be fitted to the data (e.g. the 12 quark and lepton 
masses). This doesn't mean that the theory is wrong, but suggests that there is 
probably a more complete theory out there. The electroweak theory has been ex- 
tremely successful; despite extensive experimental programs designed to test it, no 
statistically significant deviations from it have been found. 

Calculations in the electroweak theory are carried out perturbatively, in powers 
of the fine structure constant a — 1/(137.0359895 ± 0.0000061) for a momentum 
transfer of zero (at higher momentum transfers a increases somewhat; at the Z mass 



Chapter 1. Introduction 



6 



it is approximately 1/128). Because a is small, the perturbation theory converges 
reasonably well. However, for QCD, the situation is rather different. Because the 
gluons carry a colour charge, they can interact among themselves, leading to a cou- 
pling constant as that decreases for larger momentum transfer, and increases as the 
momentum transfer goes towards zero. This behaviour at large momentum transfer 
is known as asymptotic freedom, and means that perturbative QCD is valid for high 
momentum transfers, but breaks down at the lower momentum transfers where the 
coupling becomes strong. 

Unfortunately, the strong interactions found in the everyday world (the binding 
of quarks into nucleons, and of these into atomic nuclei) occur at low energies (and 
hence low momentum transfers). At present, lattice methods are the only rigorous 
way to do non-perturbative calculations with QCD. However, while great progress has 
been made using these techniques, it is not yet clear how long it will be until reliable 
results can be obtained, especially for complicated problems such as the calculation 
of decay widths of excited states. 

1.2 The Hadrons 

The six quarks and six leptons (plus their antiparticlcs) may make up matter, but 
only three of them make up the everyday matter around us. All of the matter in our 
everyday world is made up of atoms. An atom consists of a dense core, the nucleus, 
surrounded by a cloud of electrons. The structure of the electron cloud is responsible 
for the chemical properties of the atom. 

The nucleus is made up of a number of nucleons (protons (charge +1) and, usually, 
neutrons (charge 0)) bound together by the strong interaction. The atomic number of 
the atom is given by the number of protons; a neutral atom also has the same number 
of electrons. The number of neutrons is very approximately equal to the number of 
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protons. Both protons and neutrons are made up of three hght quarks: the proton 
from two m's and a d, and the neutron from a u and two d's. (Note how the charges 
of the quarks add to give the correct nucleon charges.) Everyday matter is made up 
of these three constituents (e~, u and d) because being the hghtest members of their 
'1.1), they are stable against decay^; neutrinos are also common, but because 



rows m 



they only interact via the weak interaction they interact with everyday matter very 
rarely, and we are not aware of them. 

Protons and neutrons are examples of baryons, one of the two known types of 
structures which quarks can form. Baryons are made up of three quarks {qqq)] an 
anti-baryon would be made up of three antiquarks. The other known type of structure 
is the meson, which is made up of a quark and an antiquark, qq (so an anti-meson is 
just a meson). The complicated structure of QCD means that groups of quarks can 
be bound together only for certain configurations which can have no net colour (so 
they are in a colour singlet state). It also means that the attractive force between 
coloured objects is huge, so they are always confined together into colourless objects. 
This property is called confinement. 

In a baryon each quark is a different colour, and the three colours mixed together 
produce a colourless object. In a meson the quark has a colour and the antiquark 
has the corresponding anti-colour to cancel it, again producing a colourless object. 
It may also be possible to have bound qqqqqq or qqqq states, etc., but no observed 
states have been firmly identified as such. Even more interesting would be glueballs 
(states consisting only of gluons) or hybrids (states consisting of both quarks and 
gluons - see below). Once again, no observed states have been firmly identified with 
these structures. These strongly-interacting states built up of quarks and gluons are 
collectively called hadrons. 



^Actually, a u quark can decay into a d via the weak interaction, and vice versa. However, a u 
or d always remains, so taken together they are stable. 
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In order to explain glueballs and hybrids, it helps to introduce the flux-tube 
picture of QCD. At higher energies, where perturbative QCD is valid, the picture 
of a gluon as a discrete particle exchanged by quarks is appropriate. However, at 
lower energies where the coupling becomes stronger there are indications that a more 
appropriate picture is that of tubes of chromoelectric flux connecting the quarks [0. 
These could be thought of as being similar to the lines of an electric field connecting 
electrically charged objects, but because of the self-interacting nature of the gluons 
the lines collapse down to a tube. This picture incorporates confinement - if the flux 
tubes have a certain energy per unit length, then as the quarks get further apart the 
energy of the tube would go up linearly, making it impossible to separate the quarks. 

A meson or baryon is specified by the flavours of the quarks, the space wavefunc- 
tion of the quarks, and the wavefunction describing how the quark spins combine. 
Much like a hydrogen atom, the space wavefunction can have angular and/or radial 
excitations which are identified as different meson or baryon states. It has been the- 
orized (see for example Reference 0) that the fiux-tube can also be excited, either 
vibrationally (like a vibrating string) or topologically (more complicated topologies 
of the flux-tube than just a single line). The resulting particles would be meson and 
baryon hybrids. The glueball can also be thought of in this picture as a flux-tube 
connected into a loop, with no quarks present. 

The experimental identification of hadrons other than baryons and mesons would 
be an important confirmation of our understanding of QCD. 



1.2.1 The Mesons 

As mentioned above, a meson is specified by its flavour, and its spin and space 
wavef unctions. The flavour wavefunctions that we will deal with in this work are 
given in Appendix [B.2|. The spin wavefunction is particularly simple: the two spins 



of ^ can add to give S = {a singlet) or 5* = 1 (a triplet); for more details see 
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Appendix |D.3| . Like the space wavefunction of the hydrogen atom, meson space 
wavefunctions also have a radial quantum number n and a quantum number for orbital 
angular momentum, L (but unlike the hydrogen atom, n and L are independent). The 
total angular momentum is given hj J = L + S. The space and spin state of a meson 
is often specified by spectroscopic notation, where the 25* + 1 is 1 or 3 for 

the singlet or triplet, respectively, and L is given by a letter (S, P, D, F, G, H,... being 
0, 1, 2, 3, 4, 5,...). If n = 1 it will often be neglected in the notation. It is also useful 
to give the parity and charge conjugation^ (where applicable) quantum numbers of 
the state - these are usually expressed with J as J^'^ . 



In Table |1.2| we list the known mesons constructed of the light u, d and s quarks 
(which are all that we will be concerned with in this thesis) according to the Particle 
Data Group |T]]. The rows have different space and spin states and the columns 
represent different flavours. For the 1 = 1 and I = \ column headings the quark 
states listed are considered to be different charge states of the same meson (and the 
antiparticles), because of the isospin symmetry between the u and d quarks. For the 
/ = column the two mesons in each row are listed together because the proportion 
of -^{uu + dd) and ss states in each is not the same for all of the rows. The C 
assignment in J^*^ only refers to the neutral members of the row. The names of the 
mesons are related to their flavours, the optional subscript after the name gives J, 
and the optional number in parentheses after the name is the mass of the meson in 
MeV. 

1.2.2 The Quark Model 

Because perturbative QCD does not work at low energies, and non-perturbative cal- 
culations have yet to produce detailed results, we try to calculate the properties of 



■^Parity, charge conjugation, and other conservation laws and invariance principles are discussed 



briefly in Appendix A. 4 
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jPC 


J = 1 
(-ud, -k^iuu-dd), du] 


1 = 
(-l=(uu + dd), ss) 


2 _ 

{-US, -ds, -sd, su) 




n — 1- 
U 


TT 


1], 1]' 


I/' 


1 3 c 


1 — 
i 


P 


iO, (f) 




1 1 P 


1 H — 
i 


I, / 1 oQ c;^ 




Ai(^iz /Uj, Ai(^i4UU ) ' 


1 3 p 


U ^ 




/o(1300j, /o(980j 


Ao(1430j 


1 3 p 


1 ++ 




/i(1285j, /i(1510j 


T/' /'I o7n\ /'1/lnn^t 
Ki[^lZi\]), Ai(^i4UUj' 


1 3 p 




a2(^iozU j 


/2(1270j, /2(1525j 


A2(1430j 


1 1 n 




7r2(^io/Uj 




K2[l ( (\J) 


l3n 
i i^l 


1 — 


PyL 1 UU J 


CU^^iOUU j. 


A l^iOoUJ 


13/^2 


2— 






ir2(1820) 




3^^ 


P3(1690) 


cu3(1670)i 03(1850) 


ir3*(1780) 


13^4 


4++ 


04(2040)"^ 


/4(2050), /4(2220)^ 


ir^(2045) 




0-+ 


7r(1300) 


^(1295), 


^(1460)"^ 


23^1 


1~ 


p(1450) 


cj(1420), 0(1680) 


K*(1410)^ 


23P2 


2++ 




/2(1810)^, /2(2010) 


7^2* (1980)^ 


3^So 


0-+ 


71(1770)"^ 


7^(1760)^, 


ir(1830)^ 



Table 1.2: The meson spectrum for the hght quarks. This table is mostly taken from 
the Review of Particle Properties |^ by the Particle Data Group. 
^ These states are not included in the Meson Summary Table of the Review of Par- 
ticle Properties. 

^' The Particle Data Group believes that the l^Pi and l^Pi strange mesons are 
nearly equal (45°) mixtures of the i^i(1270) and 7^1(1400). See Section and 
Appendix |A.7| . 

"I" The Particle Data Group notes that the i^'*(1410) could be replaced by the i^'*(1680) 
as the 2^5*1 strange meson. 



hadrons using models inspired by QCD rather than the full theory itself. The quark 
model^ is one such attempt, and a very successful one. 

In a quark model of a meson, the wavefunction describing the relative motion 
of the quark and antiquark is obtained by solving the Schrodinger equation with 
a Hamiltonian inspired by QCD. As an example of such a model, consider that of 
Godfrey and Isgur Their effective potential, Vqq{'p,r), contains the effects of 



^In this thesis, when we refer to the quark model, we mean the constituent quark model in 
particular. 
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a Lorentz-vector one-gluon-exchange interaction at short distances and a Lorentz- 
scalar linear interaction that models confinement. Vqq{p, r) was found by equating the 
scattering amplitude of free quarks, using a scattering kernel with the desired Dirac 
structure, with the effects of the effective potential. Because of rclativistic effects 
the potential is momentum dependent in addition to being coordinate dependent. To 
first order in (v/c)^, Vqq{'p,r) reduces to the standard nonrelativistic result: 

Vqq{p, r-) ^ V{f) = Hlf + Hlf' + Hf^ + H^^, (1.2) 

where 

Hlf = C + 6r + '^Fq ■ Fq (1.3) 

includes the spin-independent linear confinement term and Coulomb-like (from one- 
gluon-exchange) interaction. 



O fivqfiiq 



is the contact part of the colour-hyperfine interaction. 



mqlTlq 



?>Sq ■ f Sq ■ f ^ ^ 
2 '-'g ■ '-'g 



Fq ■ Fq (1.5) 



is the tensor part of the colour-hyperfine interaction, and 



SO^^SO(CM)^^SO(TP) 



is the spin-orbit interaction with 



SO(CM)^_«^ ^^^^^^^| .^^ (1.7) 



r-^ \ mqUlq mqUlq 171^ m| 
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its colour magnetic piece arising from one-gluon-exchange and 



^so(TP) ^ 1 dHgg''^ ( Sg ^ Sq . ^ 



the Thomas precession term. In these formulae, niq {rriq) and Sq (Sq) are the mass 
and spin of the quark (antiquark), L is the orbital angular momentum between the 
quark and antiquark, as{r) is the running coupling constant of QCD, and the F's are 
related to the Gell-Mann A-matrices and have expectation value {Fg ■ Fg) = —4/3 for 
a meson. 

Godfrey and Isgur's model is actually a relativized quark model, as opposed to a 
nonrelativistic quark model. They include relativistic effects by replacing the kinetic 

2 2 j 

energy term mg+mq+^ — Hamiltonian by the relativistic term Jp'^ + 



+ m? {p is the relative momentum"^ in the CM frame), by keeping m/E factors, 
and by applying relativistic smearing to the potentials. They also include the effects 
of annihilation interactions via gluons for the isoscalar mesons. 

The quark model has been highly successful at mapping the meson spectrum. 
Associated models of meson decay have also been quite successful at predicting the 
decays g |]. 

1.3 The Work of this Thesis 

We mentioned above that perturbative QCD does not work at low energies, and 
that non-perturbative calculations have yet to produce detailed results. (In fact, 
the problem arises with any strongly- interacting field theory.) Until we can make 
quantitative predictions of hadronic properties using QCD, we cannot say that we 
understand the theory, nor can we make significant headway in this important regime. 



^See Appendix A. 5 for the definitions of relative and CM coordinates. 



Chapter 1. Introduction 



13 



In order to provide some quantitative predictions for the low-energy regime, an 
industry has sprung up that calculates hadron properties by using models inspired 
by QCD, rather than the full theory itself. There are three purposes of such a pro- 
gram. The first is to shed light on the dynamics of the model, and hence (hopefully) 
on QCD, by determining which facets of the model are important in predicting re- 
sults that agree with experiment. The second is to understand the properties of the 
hadrons themselves. The third is to use these calculated properties to map the hadron 
spectrum. If we understand the meson and baryon spectra well, and there are extra 
states that do not fit into them, then these states must be something new (hybrids, 
glueballs, etc.)! Confirmation of these extra states would then feed back to the first 
purpose, understanding the dynamics of QCD. 

This thesis deals with properties of mesons calculated in the quark model, which 
is one of these attempts to approximate QCD. The meson spectrum, with its large 
number of states to be studied, provides an excellent testing ground for our under- 
standing of QCD. We examine two general areas in the quark model: models of meson 
decay, and final state interactions. 

In the first area the properties we calculate are the decay widths^ for the strong 
decays of mesons. Because meson resonances are studied experimentally by examin- 
ing the particles to which they decay, decay widths are vital in their identification. As 
an example of this we study the identity of the f 4^(2220) state - although tentatively 
identified as the ^-F4 ss meson by the Particle Data Group (see Table |1.2D , its iden- 
tity has been uncertain since its discovery in 1983. We also investigate the mixing^ 
between the -R'i(1270) and 7^^1(1400) mesons, in order to see if the dynamics of the 
model are responsible for the mixing, or whether it occurs via a separate mechanism 
that has not yet been included. This is done by comparing the models' predictions 



See Appendix A.l 
^See Appendix A. 7 



for an introduction to decay widths, resonances and cross-sections, 
for in introduction to mixing. 
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for the decay widths of these mesons to experimental data. These works have already 
been published [|], 0]. 

In the second area we investigate the effects of the final state interactions in 
the reaction 77 tttt. This is again an attempt to understand the quark model 
dynamics responsible for the interaction. In addition, if we understand the situation 
in 77 two pseudoscalar^" mesons, we can then approach the situation in 77 two 
vector mesons with more confidence; the structures seen in the cross-sections of these 
processes are not well understood and are a long-standing puzzle. 

In Chapter ^ we describe in detail the models of meson decay we are interested in, 
and the means by which we applied them. In Chapter |^ we fit the free parameter of the 
models, and investigate a number of choices that must be made in their application. 
We then apply the models to the first two chosen problems. In Chapter ^ we describe 
in detail the techniques we use for calculating the effects of final state interactions. 
In Chapter |5| we apply the techniques to our third chosen problem. In Chapter ^ we 
conclude. 



A pseudoscalar particle has = , as opposed to a scalar which has = 0+. 



Chapter 2 

Two Models of Meson Decay 



2.1 The Model of Meson Decay 

The ^Pq model p, 0, also known as the Quark-Pair Creation (QPC) model, is ap- 
plicable to OZI-allowed^ strong decays of a meson into two other mesons, which are 
expected to be the dominant decay modes of a meson if they are allowed. It can also 
be used to describe two-body OZI-allowed strong decays of other hadrons. 

In this model, meson decay occurs when a quark-antiquark pair (labelled 3 and 
4) is produced from the vacuum in a state suitable for quark rearrangement to occur, 
as in Figure |2.1| . The created pair will have the quantum numbers of the vacuum, 
J^^ = 0^"^, suggesting that they are in a ^Pq state. There is one undetermined 
parameter 7 in the model - it represents the probability that a quark-antiquark pair 
will be created from the vacuum. The pair is assumed to be created in an SU(N) 
flavour singlet, where the number of flavours N is arbitrary, since changing it will just 
rescale the value of 7 needed (see Appendix p.2|) . The rest of the model is just the 
description of the overlap of the initial meson (A) and the created pair (sometimes 



"'^See Appendix A. 6 for an explanation of the OZI rule. 
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B 



C 




B 



C 



Figure 2.1: The two possible diagrams contributing to the meson decay A —>■ BC in 
the ^Pq model. In many cases only one of these diagrams will contribute. 

referred to by 0) with the two final mesons {B,C), to calculate the probability that 
rearrangement (and hence decay) will occur. 

For the meson wavefunction we use a mock meson defined by [jlO 



2Ea E {LaMl^SaMsJJaMjJ 

X jd^pA ^nALAM^J^PA) X^aMs^ "I^A ^A 
x\qii^PA+PA)q2{^^PA-PA))- (2.1) 



Note that given our field-theory conventions of Appendix 0, the mock meson is nor- 
malized relativistically to 



{A{nA 



2Sa+1 , 



M JaMj^ 



)(P) I AiuA^'^^'LA jAMjiP')) = 2Ea 5-\P - P') 



(2.2) 



but uses nonrelativistic spinors and CM coordinates. The meson has CM momen- 
tum Pa, and the relative momentum of the qq pair pa is integrated over all values. 
The other quantities are as follows: ua is the radial quantum number; \La, M^^), 
\Sa,Ms^) and \Ja,Mj^) are the quantum numbers of the orbital angular momen- 
tum between the two quarks, the total spin of the two quarks, and the total angular 
momentum of the meson, respectively; {LaMl^SaMs^\JaMj^) is a Clebsch-Gordan 
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coefficient^; Ea is the total energy of the meson; x^SaMs^ ' ^^e appropri- 

ate factors for combining the quark spins to obtain I^a, Ms^), the flavours to obtain 
the correct meson flavour, and the colours to obtain a colour singlet, respectively; 
ipnALAMLjj>A) is the relative wavefunction of the quarks in momentum-space; mi and 
m2 are the masses of the quark and antiquark respectively; and finally, \qiijpi) ^2(^2)) 
is the basic state of a free quark and antiquark. 

Now we can consider the decay A BC . Define the S matrix 

S = I-2m 6{Ef - Ei) T (2.3) 

and then 

{f\T\i)=5%Pf-Pi) M^^'A ^'-'B ^^c (2.4) 

which gives, using relativistic phase space^ and a meson wavefunction normalized as 
in Eq. p.2| , the decay width in the CM frame 

r = 7r^i^ . ^ . \M^^A^^B^^c\\ (2.5) 

M^(2Ja + 1) Uj^Mjj„Mjc 

Here P is the magnitude of the momentum of either outgoing meson. Ma is the mass 
of meson A, S = 1/(1 + 5bc) is a statistical factor that is needed if B and C are 
identical particles, and M'^^-'a^'^-'b'^^-'c is the decay amplitude. 
For the transition operator we use 

T = -37 E (1^1-^100) /dVs d'V45'(P3+P4)3^r(^) xf-^ <i>f h\{p,) d\{p,). 

(2.6) 

Here the momenta of both the created quark (3) and antiquark (4) are integrated over 



^See Appendix A. 3 for an explanation of Clebsch-Gordan coefficients and Wigner nj symbols. 
■^Note that we have used relativistic phase space (leading to Eq. 2.5) and a relativistic normal- 



ization (Eq. pj). We will also consider other choices of phase space/normalization in Section 3.1 
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all values, with the constraint that their total momentum is zero. Other quantities 
are as follows: 7 is the one undetermined parameter in the model^; ^^(ps) and d\{p4) 
are the creation operators of the created quark and antiquark, respectively; yj^{p) = 
p'' Yi'^{9p, (j)p) is a solid harmonic that gives the momentum-space distribution of the 
created pair^; and the other quantities are as in Eq. |2.1| . Here the spins and relative 
orbital angular momentum of the created quark and antiquark are combined in a ^Pq 
state to give the pair the overall J^'^ = 0^+ quantum numbers. 

Combining Eqs. |2]l|, ^]4| and ^]6| gives for the amplitude in the CM frame (after 
doing the colour wavefunction overlap, and transforming one of the spin wavefunction 
overlaps - see Appendix 0): 



^PEaEbEc {LaMl^SaMsJJaMjJ 

X {LbMl^SbMsbIJbMj^) {LcMLaScMsJJcMj^) 
X (Im 1 -m 1 00) ixf^Ms, X^Ms^ I xIIms, xfL) 



(0B>c l0A0o ) /(P,mi,m2,m3) (2.7) 
^^|2^i4|^i2^34^ I{-P,m„m„m,) . 



The two terms in the last factor correspond to the two possible diagrams in Figure ^TT 
- in the first diagram the quark in A ends up B; in the second it ends up in C. The 
superscripts within the spin and flavour wavefunction overlaps serve to identify the 
(anti)quarks on the left side with the corresponding (anti)quarks on the right side. 



"^Our value of 7 is higher than that used by Kokoski and Isgur by a factor of a/OGtt due to 
different field theory conventions, constant factors in T, etc. The calculated values of the widths 
are, of course, unaffected. 

^The presence of the solid harmonic can be demonstrated by examining the matrix element for 
the creation of a quark-antiquark pair from the vacuum with the scalar operator constructed from 
the field: {q{pq) q[Pq) \ '^[x')'^[x") \ 0). The sohd harmonic and the Clebsch-Gordan coefficient can 
both be extracted from the spinor portion of the result. 
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The momentum space integral I {P, mi, 7712,7113) is given by 

I{P, mi, m2, ms) = J d?p rnsL.M.J,^/ + P) rncLcM.J,^/ + P) 

X^n,L,M,iP + p)yT{P) (2. 



where we have taken P = Pb = —Pc 



c- 



2.2 Calculating Meson Decay Widths with the 
Model 

In order to calculate the decay amplitudes of the ^Pq model, we use the techniques of 
Roberts and Silvestre-Brac 0. These techniques require that the radial portions of 
the meson space wavefunctions be expressible in certain functional forms, which en- 
compass simple harmonic oscillator (SHO) wavefunctions. In what follows, we assume 
that the radial portions of the wavefunctions are expressed as linear combinations of 
the first + 1 SHO radial wavefunctions (see Appendix [B.l| ), 

N 

RLAPA) = Y.daRZ=a,LAPA). (2.9) 
a=0 

The (3ai etc. found below are the oscillator parameters of the wavefunctions. 

The decay amplitudes of Eq. are converted to partial wave amplitudes by 
means of a recoupling calculation (see Appendix |^). Then the whole expression for 
the amplitudes, including the integrals of Eqs. and |E.4] , is converted into a sum 
over angular momentum quantum numbers. The graphical methods for manipulating 
angular momentum expressions found in Zare [1^ were particularly useful in this 
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regard. Our result, very similar to that of Roberts and Silvestre-Brac, is 



X E(-i)'' 

5' 



1 1 

2 2 

1 1 

2 2 

Sb Sc S' 



Sa 
1 



E 

L' 



Sb Lb Jb 
Sc Lc Jc 
S' L' S 



y^lT, ' '^^ 
L" I L" S' 1 



X < 



'sVc I0A 0o') e{LA, Lb, Lc, L, L', L", P, m,, ms, mg) + 



S' U S 
L J A L" 



(2.10) 



where the modified 9j symbol is defined in terms of the 9j symbol 
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J12 




jl 
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= <^12 -^34 "^13 "^24 J ' 
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(2.11) 



and where J = \^2J + 1, 



-F'P 



2 r>2 



e{LA, Lb, Lc, L, U , L" , P, mi, ma, m^) = 



A Pb Pc 



X E cl;- cl Ct^ {x - ^) ' {x - ^) ' {x - ly^x 
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'i2,'b,'6,'7,'8 
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Lb 




h I3 1 
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X 



'l[ + l'2 + l3 + l'^ + 2n + U+V 



(x — '"I V (x 



~r 7^0 ~i~ 



P. 



0: 



1 ^ 1 ^ 1 

Pa Pb Pc 



X = 



■P'b) p. 



PlPl + PlPh + PlPh ' 



47r(2/ + l)! 



\ (2;i + l)![2(/-;i) + l]!' 



_(-!)' / h k I 




47r 



(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 

(2.17) 



/2p(/5, /e, /7, /s, ^") = (-1)^ (2p)! k k h k E 



2(7 ^7 





2a k k 


L"+l, 



4^^(4(7 +l)(p + (7)! 

k k L" 
k h 2a 



hp+lik: k: h: k: L") = 2 {-if +\2p + l)\ k k h k E 



4^^(4(7 + 3)(p + a 



(2.18) 

1)! 



X 



V 



2(7 + 1 ^5 h 





V 



2(7+1 ^6 ^8 





^^{2p + 2a + ^)\{p 

k k L" 
k h 2(7+1 



a)\ 

(2.19) 



where 1'-^^ = Lb — h, l'2 = Lc — k, I3 = i — k, ^'4, = La — U and Dx/^^, can be extracted 
(by equating coefficients of k'^^ q'^^ {k • q)") from 



N 



\,H,v a,&,c=0 



X a/8 a! 6! c! r(a + + f ) T{h + + |) r(c + + f ) 



E E E(-i) 

ma=0 mj,=0 mc=0 



ma+nn,+mc 



PT^'PbPc - "^'^)K& - "i6)!(c - m,)! 
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^ T{ma + La + |)r(mfe + Lb + l)T{m, + Lc + |) 

ma ma-ai mt mt-bi rric rric—ci ^03+63+03 

a^o a^o to to to to a,\ a,\ a.lh.l h\ h,\ c,\ c,\ c,\ 

X (X - ^;f'^"'' (X - ^;)''^"'' ^er^^^^^^ . ^-3.3+.3-^C3^ 

(2.20) 



where 03 = nia — ai — 02, = nii, — bi — 62 and C3 = — Ci — C2. 

These expressions are generah they can be used for any meson decay where the 
radial portion of the wavefunctions can be expanded in terms of SHO radial wave- 
functions. In order to obtain symbolic results in terms of the various parameters, we 
coded these expressions into routines for the symbolic computation package Mathe- 
matica (jT^. In principle, these routines are limited only by the size of the symbolic 
problem that results, and the available computer resources - in practice, decays in- 
volving mesons with more complicated wavefunctions (i.e. having high N in Eq. p.9| , 
or high L) would have to be done numerically. The overlaps of the colour, flavour 
and spin wavefunctions of the mesons and the created pair are calculated using the 
techniques of Appendix 0. 



2.3 The Flux-tube Breaking Model of Meson De- 
cay 

In the flux-tube picture a meson consists of a quark and antiquark connected by a 
tube of chromoelectric flux, which can be treated as a vibrating string. For mesons 
the string is in its vibrational ground state. (Vibrational excitations of the string 
would correspond to a type of meson hybrid, particles whose existence have not yet 
been confirmed.) Meson decay occurs when the fiux-tube breaks at a point, and a 
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created pair (3,4) 



q (1) ' ~^ T(2) 



Figure 2.2: The position-space coordinates used in the flux-tube model. The 
cigar-shaped dashed line shows a possible surface of constant Wmm- 

quark- ant iquark pair is created from the vacuum to connect to the free ends of the 
flux-tubes, leaving a final state consisting of two mesons. 

The flux-tube breaking model of meson decay is similar to the ^Pq model, but 
extends it by considering the actual dynamics of the flux-tubes. This is done by 
including a factor representing the overlap of the flux-tube of the initial meson with 
those of the two outgoing mesons. Kokoski and Isgur |^ have calculated this factor by 
treating the flux-tubes as vibrating strings. They approximate the rather complicated 
result by replacing the undetermined parameter 7 in the ^Pq model with a function of 
the location of the created quark-antiquark pair, and a new undetermined parameter 
7o: 

7(f,w)=7oe ' (2.21) 



Here b is the string tension, where a value of 0.18 GeV^ is typically used, and ifmin 
is the shortest distance from the line segment connecting the original quark and 
antiquark to the location at which the new quark-antiquark pair is created from the 
vacuum (see Figure p.2|) : 



, w'^sin^O, ifr>w|cos6'| 

^min=<| , , • (2.22) 

r + w — 2rw I cos ^ I , if r < w I cos 9 1 
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To incorporate this into the ^Pq model, we first Fourier transform Eq. |2.8| so that 
the integral is over position-space. We then pull the parameter 7 (from Eq. |2.7| ) inside 
the integral, and replace it by the function of position 7(r, w). The expression for the 
amplitude in the fiux-tube model is then the same as that of Eq. |2.7| except that 7 is 
replaced by 70, and I {P, mi, 1712, m^) is replaced by 



/ (P, mi, ma, mg) = --— y d f d w "^Ij^LbM^I-w - r) iJjI^LcmAw - r) 

.yr{[(P + ^%.) 'P.....J:rAi^^J e-*""™" ' ('"-'^+ '"''^ (2,23) 
where m+ = — ^ 1 ^ — , m_ = — ^ — and the ih^s are now the relative 



wavefunctions in position-space. The last exponential function in Eq. |2.23| is the 



combination of the plane waves representing the CM motion of mesons B and C 
(the plane wave for A is unity, since it is at rest in the CM frame). It comes from 
Fourier-transforming the delta-functions implicit in Eq. that specify that the CM 
momenta of the mesons are held fixed to constant values. Note also that we chose to 
apply the gradient in the argument of the solid harmonic to iPuaLaMl^ only. We could 
have manipulated the Fourier-transform^ to apply it to just tpuBLsMLg "-"^ i^ncLcMLf^^ 
or some combination of all three. In fact since B and C, being lighter than A, tend 
to be less excited, choosing the wavefunctions of i? or C would in general be simpler. 



word of explanation: „ in Eq. 2.6 is rewritten as po, the relative momentum of the 



created pair. Depending on which delta functions from the anticommutators of the annihilation and 
creation operators we use to replace jjq before doing the Fourier-transform, we end up with different 
combinations of p^, ps and pc in the argument of the solid harmonic. Each pi becomes a gradient 
applied to ipuiLiML- during the Fourier-transform. 
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2.4 Calculating Meson Decay Widths with the 
Flux-tube Breaking Model 

In order to calculate the decay amplitudes of the flux-tube breaking model, there 
are a number of integrals that must be evaluated: two 3-dimensional integrals from 
Eq. p.23| and, if we choose to use the recoupling calculation to convert to partial wave 
amplitudes (see Appendix ^), another 2-dimensional integral from Eq. |E.4| . 

We wish to develop general routines to calculate any decay amplitude. This 
prevents us from making the common assumptions that can simplify Eq. |2.23| , such 
as neglecting mass differences between the quarks, restricting the orbital angular 
momenta of each meson to be or 1 so that their spherical harmonics can be expressed 
as dot products of their argument and a constant vector, or simplifying the definition 
of Wmin to rcmovc the complication that arises for r < w |cos^| (i.e. the definition 
currently used for r > w|cos0| would be used for all values). The main goal of 
simplifying Eq. p.23| is to allow more of the integrals to be done analytically, to 
reduce the computer time required for a calculation. 

In order to save ourselves two integrations, we choose to use the Jacob- Wick 
formula to convert to partial wave amplitudes (see Appendix 0). (Remember that 
this involves choosing P to lie along z, the z axis.) This leaves us with integrals over 



the variables r, 9r, 4>r, w, 9^, and 0^. The divided form of Eq. p.22| suggests that 



we will not be able to do the integrals over any of the variables on which it depends. 
This immediately rules out r and w. If r and w are expressed in the same coordinate 
system, the angle between them {6) is given by 

cos 6 = sin 6r sin 6^^ cos (0^ — 4>w) + cos 6r cos 6^ (2.24) 



which presents a problem for the other variables as well. The solution is to use the 
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freedom to choose the coordinate system of w to set Zy, to lie along f. We also choose 
to orient Xu, and i/w such that the original has (pw = in the w coordinate 

system. Then 9 = 9^], and 

f-P = rP cos 9r 

w ■ P = [sin 6^^ sin 6^^ COS 0^ + COS 6^^ COS . (2.25) 



A wrinkle with this approach is that in Eq. p. 23 we also have expressions like 



w—r 



where the vectors must be expressed in the same coordinate system. The expression 
will have to be replaced by R~^{(f)r, 9r, TT)w — f, where R{(f)r, 9r, tt) is the active rotation 
of a vector in the r coordinate system to one in the tv coordinate system. To evaluate 
the angular part of an expression such as 'iPubLbMl {~R~^{4>r,(^r,'^) w — r), we first 
use the following identity for solid harmonics 

3^r(«iri + «2r-'2) = EC,^'/^-'^ i: {hm,{l-h){m~m,)\lm)y]::^{n)ytT\^2)^ 

(2.26) 



where is defined in Eq. |2.16| . Then to handle the rotation we make use of the 
following result 

yr{R-\^r. Br. n) td) = [7^(0,, 9r, Tr)}^^!^)]* = E D^^'i^r, 9,., Tt) J^PV), 

m'=—l 

(2.27) 

where Dl^^,{(f)r, 9r, vr) is a rotation matrix. 

The final wrinkle is that the substitution for P required in the evaluation of 
/^*(— P, m2, mi, ms) requires more than just taking P —P when two different 
coordinate systems are used for f and w. The solution is to first use the symmetry 
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properties of the angular parts of the wavefunctions to obtain 



,m2,mi,m3j, 



(2.28) 



after which our particular choice for the w coordinate system can be used to evaluate 
the dot products. At this stage, the expression for the partial wave amplitudes is 
(before applying Eqs. |2.26| and p.27|) 



M^^(P)= 70 



327r(2L + 1)EaEbEc 



E 



Ml^ ,Ms^ Mhg MSb MJb' 



2Ja + 1 

x{LQS{Mj^+Mj^)\Ja{Mj^+Mj^)) {JbMj^JcMjJS{Mj^+Mj^)) 
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We can now carry out the 6r, 4>r and (pw integrations analytically, in principle. In 
practice, we will only do those for (pr and 0^. The integral over (py^ is done by means of 
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the residue theorem. The integral over (j)r depends upon the form of the wavefunction 
of meson A - we will again assume that the radial portions of the meson wavefunctions 
are expressed as linear combinations of the first + 1 SHO radial wavefunctions, as 
in Eq. p.9| . Once this assumption is made, the (pr integral does not require any special 
techniques; the final result after doing it is rather complicated but we give it here for 
completeness: 
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where (iJ„^/(^) is related to the rotation matrix by 



(2.32) 



Pi^{x) is an associated Legendre function, Ji{x) is a Bessel function of the first kind, 
and sign (a; ) is +1 if a; > and —1 if a; < 0. The g functions are given by 
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(2.35) 



where L^{x) is an associated Laguerre polynomial. 

We have automated the calculation of partial wave amplitudes with these expres- 
sions, using routines written for Mathematica [0 and FORTRAN. In the Mathemat- 
ica code, an integrand for each partial wave amplitude is prepared symbolically and 
converted to FORTRAN code. The FORTRAN program carries out the integrations 
numerically using either adaptive Monte Carlo (VEGAS [|T^) or a combination of 
adaptive Gaussian quadrature routines. Again, these routines are usable for any me- 
son decay where the radial portion of the wavefunctions can be expanded in terms of 
SHO wavefunctions, and are limited only by the size of the problem and the available 
computer resources. The overlaps of the colour, flavour and spin wavefunctions of the 
mesons and the created pair are calculated using the techniques of Appendix 0. 



Chapter 3 

Models of Meson Decay: Two 
Applications 



3.1 Setting the Parameters of the Models 

In order to use the models of meson decay introduced in Chapter |^ we must fit the 
undetermined parameter 7 (70 for the fiux-tube breaking model), and we must decide 
which meson space wavefunctions to use. In addition, we must decide which phase 
space/normalization is best suited to our purposes - we discuss this below. 

In this section we consider three types of space wavefunctions, described in some 



detail in Appendix p.l| : SHO wavefunctions with a common oscillator parameter 
(/?), wavefunctions from the relativized quark model (RQM) |^ of Godfrey and Isgur 
described in Section |1.2.2| , and SHO wavefunctions with the effective /5's of Kokoski 
and Isgur [Q. The RQM wavefunctions and SHO wavefunctions with effective /5's are 
fully determined and have no parameters left to choose. For the SHO wavefunctions 
with a common /5, we take f3 = 400 MeV, the value used by Kokoski and Isgur 
which lies safely within the range of the effective values fitted by them for the various 
mesons. We also perform another fit of both 7 and (3 simultaneously. 
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We must also decide how to normalize the mock meson wavefunctions and what 



kind of phase space to use to calculate the decay widths |T^. We investigate three 
choices. In Chapter ^ we normalized the mock meson wavefunctions relativistically 
to 2E (Eq. ^I2|) and used relativistic phase space (leading to Eq. p.5|) , which leads to 
a factor of EbEc/Ma in the final expression for a decay width in the CM frame {E 
is energy, M is mass). We will refer to this as relativistic phase space/normalization 
(RPSN). One could also argue that since our decay models are primarily nonrelativis- 
tic, it makes sense to use a nonrelativistic normalization and nonrelativistic phase 
space, which would replace the factor by MbMc/Ma- We will refer to this as non- 
relativistic phase space/normalization (NRPSN). However, there are also arguments 
p!5| that heavy quark effective theory fixes the assumptions in the mock meson pre- 
scription and suggests that the factor be replaced by MbMc/Ma, where the Mj are 
the calculated masses of the meson i in a spin- independent quark-antiquark potential 
1^.^ We will refer to this as the Kokoski-Isgur phase space/normalization (KIPSN). 

Some comments about the details of the calculations made in this chapter are 
in order. For both types of SHO wavefunctions we use quark masses in the ratio 
rriu '■ rriii : m^j = 3 : 3 : 5 - this differs from the calculations of Kokoski and Isgur 

which ignored the strange-quark mass difference. For the RQM wavefunctions 
the quark masses have already been fitted: = 220 MeV, = 220 MeV, and 
rris = 419 MeV. We treat all mesons as narrow resonances. Meson masses are taken 
from the Review of Particle Properties 1994 [Q] if the state is included in their Meson 
Summary Table^. If it is not, then the masses predicted in Ref. are used. This 
includes the masses of the 1^F2, l^F^ and l^Fi ss mesons: 2240 MeV, 2230 MeV and 
2200 MeV respectively. We ignore mass differences between members of the same 



-"^In other words Mi is given by the hyperfine averaged mass that is equal to the centre of gravity 
of the triplet and singlet masses of a meson multiplet of given orbital angular momentum L. The 
main effect of KIPSN is that the widths of decays with pseudoscalars in the final state are increased. 



The one exception was the 1 Pq ss state - see Table 3.7 
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isospin multiplet^ and average over the quoted masses, if different. Meson flavour 
wavefunctions are also taken from Ref. ||^, and are summarized in Appendix \ii.2f . 

In order to evaluate the above wavefunctions and phase space/normalization 
schemes and to fit the value of 7 (70) we carried out a series of least squares fits 
of the model predictions to the decay widths of 28 of the best known^ meson decays; 
we minimized the quantity defined by = J2i{^T°'^'^'' ~ '^T^Y / '^Vi "where ar, is the 
experimental error. ^ The 28 decays used, and the experimental values of their widths 
, are shown in Table |3.1] . 

The ratios of the model predictions to the experimental values for the widths of 
the 28 decays are plotted in Figures ^7\\ to |3.5| , on logarithmic scales. The fitted values 
of 7 (70) and the resulting per degree of freedom (dof) are listed in Table [3^. We 
do not quote errors on the fitted values of 7 (70) because they are negligible (being at 
most a change in the last digit quoted) compared to the inherent uncertainties of the 
models. We have performed fits for the ^Pq model with both types of SHO wavefunc- 
tions, and for the flux-tube breaking model with all three types of wavefunctions. We 
did not use the RQM wavefunctions with the ^Pq model because the large number (up 
to 13) of radial SHO wavefunctions in the linear combination (see Eq. p.9| ) exceeds 
the capacity of our symbolic routines on the available computer hardware^. Because 
the results of the ^Pq and flux-tube breaking models appear to be very similar, we 



''The one exception was for the decay </) K^K^ where the charged and neutral kaon mass 
difference is significant to the phase space. 

"^For the 77 and 77' we assumed perfect mixing, corresponding to 9p — —9.7°, but we also considered 



the effect that using Op = —20° would have on the results of Tables 3.5 through 3.9. We found that 
while individual decay widths varied by up to a factor of approximately 2, the changes to the total 
widths were negligible compared to the other uncertainties of the calculations. 

'"^We also required that the decays we selected not involve any particles that undergo significant 
mixing (e.g. 77, 77') in order to minimize uncertainties in the calculation. 

^For the calculations in the flux-tube breaking model, a 1% error due to the numerical integration 
was added in quadrature with the experimental error. 

^This problem arises for the ^Po model routines, which are symbolic (and hence are slower and 
require more memory), but not for the flux-tube breaking model routines, where the brunt of the 
work is done numerically. 
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Table 3.1: The 28 meson decays used in our fits of the decay models' parameters. 
The experimental widths are shown, as are the model results for six combinations of 
decay model, wavefunctions, and phase space/normalization (specified in the second, 
third and fourth rows of the heading, respectively) . 
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Figure 3.1: Ratios of the predictions of the ^Pq model of meson decay to the exper- 
imental values, for the widths of the 28 decays used in the fits. In this figure SHO 
wavefunctions with a common (3 were used. The three graphs use different phase 
space/normalization schemes. 

chose not to write a numerical version of the ^Pq routines to enable us to use the 
RQM wavefunctions. For each decay model and wavefunction combination, we tried 
all three phase space/normalization schemes. 



When examining Figures |3.1| to p.5|, note that if a particular model described the 
data exactly, all the points for it would lie along the vertical solid line at rm°deYpexp _ -^^ 
Because 7 (70) is an overall factor in the amplitude, changing it would move all of the 
points to the left or right by the same amount. Some distributions appear lopsided 
about rmodey-pexp _ because the experimentally well known decays dominate the 
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Figure 3.2: Ratios of the predictions of the ^Pq model of meson decay to the exper- 
imental values, for the widths of the 28 decays used in the fits. In this figure SHO 
wavefunctions with effective /9's were used. The three graphs use different phase 
space/normalization schemes. 

fit. 



The first thing to notice from the figures and the values of Table |3.2| is that the 
results with NRPSN are significantly worse than the ones using the other two choices 
of phase space/normalization. We will hence not use NRPSN for our calculations. 
Similarly, the SHO wavefunctions with effective /?'s do worse than either of the other 
two types of wavefunctions, and we will not use them either. That leaves us with six 
different combinations of decay model, wavefunction, and phase space/normalization 



these are the combinations given in Table 3.1. Some general comments are in order. 
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Figure 3.3: Ratios of the predictions of the flux-tube breaking model of meson decay 
to the experimental values, for the widths of the 28 decays used in the fits. In this 
figure SHO wavefunctions with a common f] were used. The three graphs use different 
phase space/normalization schemes. 



First, the values of Table p.2| make it clear that these models are not very 
accurate^ - the best they can hope for is to predict a decay width to within a factor 
of 2, and even larger deviations are common. Since they are coarse models of a 
complicated theory, this is not surprising. In addition, the six combinations can give 
rather different results. We will use this to our advantage by taking the spread in the 
predictions of the different combinations as a rough guide to the trustworthiness of 



^It should be noted that past practice has often been to look for agreement in the amplitudes, 
rather than the widths, which tends to make the situation look better. 
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Figure 3.4: Ratios of the predictions of the flux-tube breaking model of meson decay 
to the experimental values, for the widths of the 28 decays used in the fits. In this 
figure SHO wavefunctions with effective /?'s were used. The three graphs use different 
phase space/normalization schemes. 

the calculations. 



From Table |3.1| and Figures |3.1| through |3]^ one can see that the results for the 
^Po and flux-tube breaking models for the SHO wavefunctions are very similar, as 
previously noted. ^ This is in agreement with the observations of Kokoski and Isgur 
0, who argued that the overlap of the space wavefunctions in the ^Pq model makes 
it unlikely that the decay will proceed if the quark-antiquark pair is created far from 



^The one exception to this is the S-wave decay i^Q (1430) Ktt which seems particularly sensitive 
to the model. 



Chapter 3. Models of Meson Decay: Two Applications 



39 



RPSN NRPSN KIPSN 



bi(1 235) ^cjTT 
0^(1 320) ^pTT 

0^(1 320) ^KK 
7x2(1 670) ^f2(1270)TT 
7x2(1 670) ^pn 
7x2(1 670) K'(892)K 

P3(1 690)^7TTT 
P3(1 690) ^WTT 

P3(l 690) KK 

(2(1270) ^TTTT 

(2(1 270) ^KK 
f^(2050)^ww 

(4(2050) ^7T7T 

(4(2050) ^KK 

(p K*K" 

(2^(1 525) ^KK 
K'(892) K7\ 

Ko"C430) ^Ktt 

Kj'C -150) Ktt 
K2X143O) K'(892)7T 

K2"(1430)^Kp 

K2'(1430) Koj 

K3'(1780)^Kp 
K3XI 780) K"(892)7T 

K3'(1 780) Ktt 

K4'(204b) Ktt 
K4X2045) K'(892)99 

10"^ 10"^ 1 10 10"^ 10"' 1 10 10"^ 10"' 1 10 

pmadel / pexp 



Figure 3.5: Ratios of the predictions of the flux-tube breaking model of meson decay 
to the experimental values, for the widths of the 28 decays used in the fits. In 
this figure RQM wavefunctions were used. The three graphs use different phase 
space/normalization schemes. 



a straight line connecting the two mesons. This would have the same effect as the 
functional form of 7(r, w) in the ffux-tube breaking model. Despite the similarity in 
the individual decays, the values of the fits are better for the fiux-tube breaking 
model, indicating at least some preference for this model. 

In comparing the SHO wavefunctions (we will henceforth only speak of the SHO 
wavefunctions with a common (3) to the RQM wavefunctions, we note that the RQM 
wavefunctions do better for RPSN, but worse for KIPSN, showing no clear preference. 

KIPSN gives better overall fits to the data than RPSN. However, certain decays. 
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SHO (common f3) 
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16.5 
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371 MeV 
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Table 3.2: The details of the various fits of the decay models' parameters, and the 
resulting x^/dof. For the last three fits, both 7 and (3 were fit simultaneously. 

/s:^(1780) Kp and /4(2050) lulu for example, are fit better using RPSN. Also, if 
we wished to remove the dominant effect that the experimentally well known decays 
have on the fit, we could follow Geiger and Swanson fl^ and look at the spread of the 
values of 7 needed to match the model predictions to the experimental values exactly; 
if 7 is the mean of these 7 values, and the standard deviation of the distribution, 
then cT^/7 provides a measure of the spread of the results that ignores the experimental 
errors. Using this measure, RPSN is consistently better than KIPSN, as can be seen 



in Table 3.3 



We can make more specific observations using the SHO wavefunctions: Decays 
with two pseudoscalars in the final state tend to do better with KIPSN, but KIPSN 
generally underestimates decays of high L mesons with vector mesons in the final 
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27.6 
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KIPSN 




17.1 
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Table 3.3: The values of 7, and for the 28 decays used in the fits, where the 7 
(70) value for a decay matches the model prediction to the experimental value. Thus 
cr^/7 provides a measure of the spread of the results that ignores the experimental 
errors. 



state. On the other hand RPSN tends to overestimate decays with two pseudoscalars 
in the final state. Trends can also be observed with the flux-tube breaking model 
using the RQM wavef unctions. Having said all this we stress that these are only 
general observations and exceptions can be found to any of them in Table p.l| . One 
must therefore be very careful not to take any predictions of these models at face 
value but should try if possible to compare the predicted decay to a similar one that 
is experimentally well known. 

Finally, we consider the sensitivity of our results to which we must choose 
for the SHO wavef unctions. Different quark models predict different values of /?, 
so there is no clear theoretical prediction to draw upon. It should be noted that 
although meson decay widths can depend strongly on the value of /3, to some extent 
we alleviate this effect because we fit 7 to experimental data after selecting a value 
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Figure 3.6: Ratios of the predictions of the ^Pq model of meson decay to the exper- 
imental values, for the widths of the 28 decays used in the fits. In this figure SHO 
wavefunctions with a common j3 were used, and both P and 7 were fit simultaneously. 
The three graphs use different phase space/normalization schemes. 

of (3. In the fits discussed above, we took j3 = 400 MeV which is the value used by 
Kokoski and Isgur 0. In addition to these fits, we performed simultaneous fits of 
both 7 and /? to the 28 decay widths - see Figure |3]^. The overall fits improved, but 
as usual some widths were in better agreement and some in worse agreement with 
experiment when compared to the fits for (3 = 400 MeV. In general we found that the 
deviation in the results lies within the overall uncertainty we assign to the models, 
and so we choose to follow the literature and use (3 = 400 MeV. 



Chapter 3. Models of Meson Decay: Two Applications 



43 



3.2 The /4(2220) Reexamined - What Is it Really? 

The /4(2220) (formerly known as the ^(2220)) state was discovered by the Mark III 
Collaboration pf in 1983, in the decay e+e" J/tp f 4(2220) -fK+R- - 
see Figure It attracted considerable attention because of its narrow width of 
approximately 30 MeV, which was thought to be too small for a conventional meson 
of that mass. Over the years this has led to speculation that it might be a Higgs boson 



17|, a bound state of coloured scalars [Wl, a four quark state M, BO], a AA bound 



state 1 21], or a meson hybrid or a glueball [^, among others (we only reference early 
works on each suggestion). 

The /4(2220) has since been seen by other experiments, the results of which are 



summarized in Table |3.4| . The rather broad state seen by DM2 is unlikely to be the 
same state seen by the other experiments. The recent measurements by the BES 
Collaboration are of particular interest, because the previously unseen tt+tt^ and pp 
decay modes indicate a decay pattern that is roughly flavour-symmetric, in keeping 
with the glueball interpretation (see for example, the recent paper by Chao p3|)- "^^^ 
/4(2220) has not been seen by pp experiments, but the limits set are not stringent 
enough to exclude the meson identification (for example see Barnes et al. ]]23]). 



The intermediate state J / ip and the final observed states indicate that the parity 
and charge conjugation quantum numbers of the /4(2220) are and that its spin J 
must be even, although as can been seen from Table |3.4| the actual value of J is rather 
uncertain. The Particle Data Group ]jl| has tentatively identified the /4(2220) as the 
^F^ ss meson (hence its presence in Table |1.2|) , but notes that this needs confirmation. 

In order to determine whether the /4(2220) is a high-spin meson or something new 
and much more interesting, we need theoretical predictions of the properties of the 
various options in order to make the identification. In this work we carefully examine 
the decays of the and ^^4 ss mesons in order to see if the meson identification is 
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Figure 3.7: The original data from the discovery of the /4(2220) by the 
MARK III Collaboration |16|. The distribution of the K^K~ invariant mass in 
e^e~ — s> J/-?/^ ^ 7 K^K~ is shown, revealing three resonances. The errors shown are 
purely statistical. 

tenable. These are the best meson candidates because of their J^'^^s of 2"*"+ and 4+"*" 
respectively, and because their masses have been predicted to lie in the right range 
(Godfrey and Isgur 0] calculated masses of 2240 MeV and 2200 MeV respectively). 

The decay widths of these mesons have been calculated previously. The first was 
by Godfrey, Kokoski and Isgur |^| using the fiux-tube breaking model (and somewhat 
different masses: 2270 MeV and 2210 MeV respectively). Despite the fact that one 
naively expects heavy mesons to be broad, since more phase space and decay channels 
lead to larger widths, they found the total widths to be approximately 60 MeV and 
50 MeV respectively. Based on production rate arguments, they concluded that the 
/4(2220) was most likely to be the ss state. However, the analysis was not 
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Experiment 


Process 


Mass 
(MeV) 


Width 

(MeV) 


J 


SLAC-SP-032 
(Mark III) ^ 


e+e- J/ij ^7/4(2220) 


2230 ± 6 ± 14 
2232 ± 7 ± 7 


26+?« ± 17 

istll ± 10 




(GAMS-2000, 
GAMS-4000) ^ 


TT^p /4(2220) n 


2220 ± 10 




> 2 


SERP-E-147 
(MIS) ^ 


Tj-p f 4(2220) n 


2230 ± 20 


80 ±30 


2 


SLAC-E-135 
(LASS) ^ 


K-p /4(2220) A 
K+K- 


2209^1^ ± 10 


60ir 


4 


RPP Average ° 




2225 ± 6 


38+}^ 


4? 


(DM2) ' 


e+e~ ^ J/-?/; — > 7? 


2197 ± 17 


201 ±51 




(BES) s 


e+e- ^ J/V' ^7/4(2220) 

^ K+K' 

c_i. pp 
^ TT + TT- 


22301? ± 12 
22321? ± 15 
2235 ± 4 ± 5 
2235 ± 4 ± 6 


20tlt ± 12 

20^11 ± 10 
151^2 ^ 9 

19l|? ± 12 





Table 3.4: Summary of measurements of the /4(2220). The RPP average is for the 
first four entries in the table only; they feel that the /4(2220) still needs confirmation. 
In the case of the DM2 results, the resonance seen was not identified with the /4(2220), 
and so it is not included in the RPP average; the BES data is too recent. 



^ Reference JIB 
^ Reference 
^ Reference [|6 

Reference 
® Reference [||]. 
^ Reference 



^ Reference 29 . 



exhaustive in that it did not calculate the partial widths to all possible final states. 
In particular it made the assumption, which we will see to be incorrect, that the 
decays to an L = 1 meson and a K or rj were small on the basis of phase space 
arguments alone. 



•"in addition, they note that the decay to K* (892) K* (892) is dominated by the 1, D partial wave 
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In addition, Pakvasa, Suzuki and Tuan [^] have used Regge theory to estimate 
the partial width of the ss state decaying to two ^5*0 mesons as approximately 
40 MeV, with the total width expected to be a few times larger. And finally, Ono, 
Pene and Schoberl have used the ^Pq model to calculate the partial width of the 

ss state (assuming a mass of 2200 MeV) decaying to ^5*0 and ^Si states to be 
417 MeV (with the width to K* (892) K* (892) dominant at 348 MeV). In light of the 
new data, and these contradictory theoretical results, a new, complete calculation is 
required. 

We have calculated the and ss meson decay widths using the six combina- 



tions of model, wavefunction, and phase space/normalization discussed in Section 3A 



In order to judge the reliability of the results, it is useful to examine how well our 
calculations do for similar mesons. We do not have any good candidates for the ^^2, 
but the /4(2050) and ^^4(2045) differ only in flavour from the ^^4 ss. Examining Fig- 
ures |3.3| and we see that for the five decays in question, the choice of phase 
space/normalization makes a significant difference: RPSN tends to overestimate the 
partial widths, while KIPSN tends to underestimate them. For the K^{1780), which 
one would expect to behave similarly to the ^-F4 states (since both are ^L^+i), the 
situation is unfortunately less clear. We can also examine the total width of the 
/4(2050) and ^^4(2045) states. In Tables |3]^ and |3.6| we give the partial widths for all 



significant, kinematically-allowed decays and their sums, in the ^Pq model and with 
SHO wavefunctions only. It can be seen that once again the RPSN results are high, 
and the KIPSN results low. The main conclusion we can draw from these results 
is that the widths of the ^^4 ss meson probably lie between the RPSN and KIPSN 
estimates but it is difficult to guess if they are closer to the lower or upper value. 

we 



Turning now to our results for the F2 and F4 ss mesons, in Tables |3.7| and ^ 



(using the 5, L notation of Table ^^). Unfortunately, S = 1 is forbidden by conservation of parity 
and charge conjugation, making the result for this decay questionable. 
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Decay 
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0.4 
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-'-•'^-0.4 


5.4 


1.6 
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2.7 


0.8 


[KKi{1270) + c.c 






2.3 


1.2 


[K* (892) K* (892)] 


2,-D 




7.3 


2.1 


Total / J2i Ti 


208 ± 13 


435 


109 



Table 3.5: The calculated partial decay widths and total width of the /4(2050) meson. 
The available experimental results ^ are also shown. The model results are given 
for two combinations of decay model, wavef unctions, and phase space/normalization 
(specified in the second, third and fourth rows of the heading, respectively). We 
have calculated the partial widths of all kinematically-allowed, OZI-allowed, two-body 
strong decays, but only show those partial widths that are > 1 MeV in at least one 
model. For this reason the total widths may not equal the sum of the partial widths 
shown. The subscripts on the decays refer to the S and L (see Appendix ^ of the 
given partial wave - the L is in spectroscopic notation {S, P, D, F,G, H). We have 
used the meson mixings given in Appendix [A.7| where applicable. 
^ This number is the total for all partial waves. 



give the partial widths for all significant, kinematically-allowed decays and their sums, 
for all six combinations of decay model, wavefunction, and phase space/normalization. 
For a degree of completeness, in Table |3^ we also give the results for the other L = 3, 



spin-triplet ss meson, the ^F^, in the ^Pq model and with SHO wavefunctions only. 
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Total / EiTi 
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Table 3.6: The calculated partial decay widths and total width of the K% (2045) meson. 
The available experimental results [Q are also shown. For additional comments, see 



Table 3.5 



This number is actually for the final state ir*(892)7r7r, and is the total for all partial 
waves. 

For the ss state the main decay modes are, in descending order, including 
charge conjugate pairs: 



ii'i(1270)ir, fs:i(1270)ir*(892), K*{892)K, K* (892) K* (892), 

k;{u30)k, kk, /i (1510)77. 
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Decay 

3 p „^ , 
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0.4 
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5.9 
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99 
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0.7 
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0.9 
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0.3 
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0.3 
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0.3 




1046 


391 


1058 


406 


2181 


797 



Table 3.7: The calculated partial decay widths and total width of the ^^2 ss meson. 
We do not include a decay to /o(980)/o(980) because we question its assignment as 
a ^Pq meson. At a more likely mass for the ^Pq ss meson, this decay is kinematically 
inaccessible. For additional comments, see Table |3.5| . 
^ Kr is our notation for the first radial excitation (2^5*0) of the K. 
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44 


184 


46 
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72 
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2.0 
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15 
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0.6 
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20 
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21 
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31 


9.5 
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Table 3.8: The calculated partial decay widths and total width of the ss meson. 



For additional comments, see Tables |3l5| and BT? 



It is clearly a mistake to neglect the decays to an L = 1 meson and a or 77; by 
far the largest partial width is for the decay to -ft'i(1270)-S', and many of the other 
decays involving L = 1 mesons are significant. The results using RQM wavefunctions 
are surprisingly large - about double those using SHO wavefunctions. Turning to the 
results with SHO wavefunctions, we see that if we believe that the KIPSN results give 
a lower bound for the as well as the ^^4, we can expect a total width >, 400 MeV, 
in keeping with naive expectations. Even if this width is too large by a factor of 2, 
the /4(2220) cannot be the ss meson. Indeed, the ^^2 ss may well be too broad 
to be found experimentally. 

For the ^^4 ss state the main decay modes are, in descending order, including 
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Decay 

3^ SS — > 


r (MeV) from Models 


3p 


SHO 


RPSN 


KTPSN 


[A [QyzjK + c.c.\i^o 


1 /I 


oo 


[K [p't)Z)K + C.C.Ji^G 


yo 


z4 


[A (^i4iUjA + C.C.Ji 


z.o 


i.O 


[A (^i4iUjA + C.C.Ji^G 


2.2 


1.1 


[AQ(^i4oUJA + C.C.Jo,F 


5.3 


2.4 




71 


30 


[fC*(1430)ir + c.c.]2',P 


305 


137 


[K*{892)K*{892)]2,D 


66 


16 


[K*l892)K*l892)]2,G 


27 


6.6 


[Ki{1270)K*{892) + c.c.]2,P 


157 


63 


[/^(1525)r/]2,p 


73 


37 




10 


3.2 




974 


360 



Table 3.9: The calculated partial decay widths and total width of the ss meson. 
We do not include decays to /o(980)?7 or fo{980)ri' because we question the assignment 
of the /o(980) as a ^Pq meson. At a more likely mass for the ^Pq ss meson, this decay 



is kinematically inaccessible. For additional comments, see Table pT5 
charge conjugate pairs: 



ir*(892)ir*(892), KK, K*{892)K, 



Note that the total width of the ss is quite a bit less than that of the ss . This 
can be understood in terms of the orbital angular momentum between the mesons in 
the final state. Because of the higher spin of the ^^4, it tends to yield final states with 
higher orbital angular momentum in order to conserve angular momentum. These 
decays are then suppressed by the higher angular momentum barrier. Examining 
Table |3.7| , we see that the lowest orbital angular momentum final states in ss 
decay are in P-waves. All of these decays are relatively broad but the decay to 
i^ri(1270)^ is the P-wave decay with the largest available phase space. In fact, one 
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could almost order the P-wave decays using phase space alone. The analogous decay of 
the ss is in an F-wave and therefore is suppressed by a larger angular momentum 
barrier. The lowest angular momentum partial wave for ^F^ ss decays is a D-wave 
which, although it has the largest partial width of all ss decays, is still smaller 
than the P-wave ss decay. 

Looking at the total width of the ss state, if we accept that the true result 
lies between those obtained from RPSN and KIPSN, then we would expect the width 
to be 140 MeV and < 600 MeV. While this is not very specific, we are really 
interested in the lower bound, and whether this could be the /4(2220). It is possible 
to look at the behaviour of the similar /4(2050) and ^^4(2045) (and perhaps K^{1780)) 
decays in selected model, wavefunction and phase space/normalization combinations. 
For example, in the ^Pq model with SHO wavefunctions and RPSN, the decays to two 
^Si mesons are well predicted, indicating that the largest partial width of the ^^4 ss 
(to i^*(892)^*(892)) may well be close to the calculated value. On the other hand, 
the partial widths to two ^Sq mesons may be less than predicted, and to a ^So and a 
^5"! slightly more. Based on this, we think it possible, but unlikely, that the calculated 
^F^ ss width is high by a factor of 2. Indeed, we would expect the actual width to be 
higher, i.e. > 140 MeV. Unfortunately, the uncertainties associated with the models 
we are using make it impossible to rule out the identification of the /4(2220) as the 
^F4 ss meson with any surety; however, we feel it is unlikely. 

Aside from the total width of the ^F^ ss, the roughly flavour symmetric nature of 
the /4(2220) decays observed at BES suggest that it is not an ss meson, since produc- 
tion of TTTT and pp would be OZI-suppressed. We thus propose a second explanation 
of what is being observed in this mass region - that two different hadron states are 
present: a narrow state produced in J/ip radiative decay and a broader state seen 
in hadron beam experiments. The broader state would be identified with the ^^4 ss 
state. The measured width is consistent with our predictions and the LASS CoUabo- 
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ration ||2^ shows evidence that its quantum numbers are J^^ = 4"*"+. We would then 
identify the narrow hadron state observed in the gluon rich J/ip radiative decays as 
a glueball |2^. The narrow state is not seen in hadron beam production because it 
is narrow, is produced weakly in these experiments through intermediate gluons, and 
is hidden by the ss state. Conversely, the broader state is not seen in J/ip radiative 
decays because this mode preferentially produces states with a high glue content. 
Crucial to this explanation is the experimental verification of the BES results on the 
flavour symmetric couplings of the state produced in J/ip radiative decay and the 
observation of other decay modes for the broader state (primarily /^* (892)^* (892) 
and K*{892)K) in addition to the theoretical verification that the predicted tensor 
glueball is as narrow as the observed width. 

We conclude that the /4(2220) cannot be the ss meson, and that taking all of 
the data into account, it is unlikely to be the ss meson alone. Instead we suggest 
that the broad state seen in hadron beam experiments is the ss meson, and the 
narrow state seen in J/ ip radiative decay is a glueball. Further experimental data is 
needed to finally identify the /4(2220). 



3.3 Constraining the Ki Mixing Angle 



As discussed in Appendix |A.7| , the ^Pi and ^Pi strange mesons mix to produce the 



physical states i^'i(1270) and i^'i(1400). A possible source of the mixing can be seen 
in the spin-orbit portion of the Hamiltonian of Godfrey and Isgur's relativized quark 
model H], described in Section |1.2.2|. We can rewrite Eq. Ol as 



(3.1) 
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where 



H 



S0+ 



99 



H. 



so- 

99 



— + — F,-F, 



1 ai/g^f ^ 
4r 9r 



conf ^ 



(5. 



S,) ■ L. 



(S, + S,) ■ L, 

(3.2) 



Using this quark model Hamiltonian, the mass formulae for the P-wave mesons 



are 



MCPo 



1 

To' 



Mo + \{H^r) + UKT) - (^9^^) 



99 



99 



Mo + 7(^r ) - (HIT) - 2(^aT^)' 



99 



99 



M,-\{Hlf') 



(3.3) 



where the (H^g) are the expectation values of the spatial parts of the various terms 
of the Hamiltonian, Mq is the center of mass of the multiplet, and we have adopted 
a phase convention corresponding to the order of coupling L x Sg x Sg. It can be 
seen that Hg^~ leads to mixing between the ^Pi and ^Pi states. However, H^^~ is 
non-zero only if rriq ^ m^, so the strange mesons are the lightest states where this 
mixing can occur. 



If we use Eq. A. 5 to describe the mixing, we find that to diagonalize the ^Pi-^Pi 
Hamiltonian we must have 



tan(2^ 



(3.4) 



where 6k is the Ki mixing angle. In this particular quark model, the expectation 
values of the spatial parts of the various terms for strange mesons are Mq = 1378 MeV, 
= 33 MeV, (HlT) = 56 MeV, {H^P+) = 47 MeV and {H^P-) = -1 MeV 
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[ 34 . This gives rise to a prediction of 9k = 

In this section we compare the predictions of the ^Pq and flux-tube breaking 
models of meson decay to experimental data for five decay widths and two ratios of 
D to S amplitudes'^, in order to determine 6k- This mixing angle is of interest both 
for its own sake, and for what it can can tell us about the quark model Hamiltonian. 

The i\:i(1270) and Ki{UOO) can both decay to Kp and if*(892)7r; the i\:i(1400) 
can also decay to Ku. All of the decays occur in the partial waves 1, S and 1,D {in S,L 
notation). Although decays to other final states are observed they lie below threshold 
and proceed through the tails of the Breit-Wigner resonances making the calculations 
less reliable. In addition to these five decay widths, there are also experimental results 
available for the ratio of D to S amplitudes for both i^i(1270) and i^i(1400) decaying 
to K*{892)n. 

We have calculated these seven quantities using the six combinations of model. 



wavefunction, and phase space/normalization selected in Section |3Tl| . The Ki(1270) 



quantities are plotted as a function of the Ki mixing angle in Figure \i.S[ and the 



i^'i(1400) quantities are plotted in Figure |3]^ - the experimental results are also 
shown in both figures. 

In order to determine the best fit values of 6k, we have carried out a least 
squares fit for each of the six combinations of decay model, wavefunction, and phase 
space/normalization. The best fit values of 6k and their corresponding x^/dof are 
shown in Table |3.1CI| , as are the resulting values of each of the seven quantities being 



considered, and their experimental values. Examining the results, we see that as 
usual, the values of x^/dof are high, but not as high as were obtained in the 28-decay 
fits of Section |3.1| . This is not unexpected, since the similar decay 61 (1235) ujtx was 
described reasonably well in those fits. It is also clear that the SHO wavefunctions 



11 By "ratio of D to S amplitudes" we mean the ratio of the L = 2 partial wave amplitude for a 
decay to the L = partial wave amplitude for the same decay. 
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Figure 3.8: Graphs of the i^'i(1270) partial decay widths and ratio of D to S am- 
phtudes considered in this work, vs. the mixing angle 6^. The shaded regions 
are the l-a limits on the experimental measurements of the quantities In 
the key the labels refer to the meson decay model, the wavefunctions used, and 
the type of phase space/normalization used. Only three curves are shown for 
\AD/AsiKi{1270) K*(892)7r)|, because the differences between RPSN and KIPSN 
cancel in the ratio of amplitudes, making the RPSN and KIPSN curves for the same 
model and wavefunctions identical. 

do better than the RQM, and in terms of phase space/normalization, KIPSN does 
better than the RPSN. 





Value 


Value Calculated from Models of Meson Decay 


Quantity 


from 






Flux-tube 


Breaking 






Experiment 


SHO 


SHO 


RQM 






RPSN 


KIPSN 


RPSN 


KIPSN 


RPSN 


KIPSN 


Fitted mixing angle Ok 




48.4 ± 1.1° 


45.2;i:f 


46.8;i:i° 


44.2 ± 1.3° 


55.4i[ii° 


50.6^|:f 


xVdof 




38.9 


5.3 


26.0 


3.3 


231 


51.2 


r(iri(1270) ^ Kp) 


38± 10 


127 


76 


115 


70 


213 


121 


r(ii'i(1270) is:*(892)7r) 


14.4 ±5.5 


24 


13 


20 


12 


76 


34 


\AD/As{K^{m{)) ^ K*{m2)'K)\ 


1.0 ±0.7 


0.65 


0.86 


0.77 


0.99 


0.31 


0.41 


V{Ki{lAm) Kp) 


5.2 ±5.2 


15 


12 


17 


12 


10 


7.0 


T{Ki{Um) K*{m2)'K) 


164 ± 16 


350 


221 


306 


197 


653 


402 


\Aj,/As{K^{imQ) ^ K*{m2)'K)\ 


0.04 ±0.01 


0.071 


0.054 


0.069 


0.053 


0.080 


0.061 


r(iri(1400) ^ Kuj) 


1.7± 1.7 


4.3 


3.5 


5.0 


3.8 


2.7 


2.0 



Table 3.10: The partial decay widths and ratios of D to S amplitudes used in our fits of the i^i(1270) and ^1(1400) 
mixing angle 9k- The experimental values [1] are shown, as are the model results for the six combinations of decay 
model, wavef unctions, and phase space/normalization (specified in the second, third and fourth rows of the heading, 
respectively). For each combination we also give the fitted value of Ok-, and the corresponding x^/dof. The partial 
widths are given in units of MeV. All errors shown are l-o". 
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Figure 3.9: Graphs of the i^'i(1400) partial decay widths and ratio of D to S am- 
phtudes considered in this work, vs. the mixing angle 6k- See Figure for the 
explanation of the line labels and additional notes. 

The resulting values of 6k are all fairly similar, falling in the range 44 — 56°. 
The lowest values of x^/dof occur for fits near 44°, and so we take that as our 
approximate result. What can this tell us about the quark model Hamiltonian? A 
mixing angle of 44° implies {H^P~) = —142 MeV, rather different from the quark 



model value. However, the value of {H^^ ) is very sensitive to the angle in this region 
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(e.g. 0K = 41° implies {H^f') = -35 MeV, while 0k = 45° implies {H^f') = -oo), 
so the only conclusion we can really draw is that the quark model value is very low, 
if {H^^^) IS responsible for all of the mixing. We cannot rule out the possibility 
that some other interaction is partly responsible for the mixing, such as mixing via 
common decay channels. 

To consider the SO Hamiltonian further, define m = |("^g + "^g) ^-nd A = 
(mq — mq). Next expand l/rUq — l/m+A/(2m^) + ... and l/rUq — 1/m — A/(2m^) + ... 
to first order in A/m (each of these truncated results has less than an 11% error for 
the quark masses we used). In terms of these, the contributions to the spin-orbit 
Hamiltonian become 

The {{Sq ± Sg) ■ L) factors are roughly the same size, and the A/m term is not 
too far from 1 (it is > | for the quark masses we used), so the factors in the large 
parentheses contribute at the same order of magnitude. Note that the contributions 
from the Thomas precession term are equal, but those from the colour magnetic term 
differ by a factor of 2. The model gives a large {H^^~^) of 47 MeV, but a small 
(H^^^) of —1 MeV arising from a delicate cancellation of the colour magnetic and 
Thomas precession terms. It is thus conceivable that a small change in the model 
could lead to a substantially larger {H^?'~), so we are unable to draw any conclusions 
about whether the mixing is largely due to another mechanism, or the quark model 
Hamiltonian just needs some retuning. 



Chapter 4 

Effects of Final State Interactions 



4.1 Introduction to Final State Interactions 

Final state interactions (FSI's) refer to additional interactions between particles in the 
final state that are not included in the basic interaction.^ For example, in Chapter ^ 
we consider the basic interaction 77 tt+tt^, which proceeds through a scalar QED 
vertex. However, there will also be a QCD interaction between the final state pions. 
In Chapter |^, we attempt to calculate the effects of the FSI for this process. In this 



chapter, we introduce the techniques that we will apply to the problem of FSI's [p3 |. 

The techniques involve nonrelativistic quantum mechanics (including the Schro- 
dinger equation), corrected to include relativistic phase space where possible. The 
basic interaction may still be calculated relativistically: e.g. 77 tttt. This mixture 
of nonrelativistic and relativistic elements is typical of the quark model. 



^The techniques introduced here could also be used to describe additional interactions between 
particles in the initial state. However, since most physically interesting problems involve FSI's, the 
term "final state interactions" has become the name of choice. 
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4.2 Solving the Schrodinger Equation with the 
FSI Potential 

Our first step in calculating the effects of an FSI is to find the wavefunction describing 
the state of the final state particles due to the presence of the FSI only. We assume 
that the FSI can be described by a potential between the final state particles, and 
that the wavefunction can be found by solving the Schrodinger equation. 

The relative wavefunction for two bodies in a potential that only depends on 
their relative position f = fi — f2 can be found from the Schrodinger equation for a 
central-force potential 

V^. - 2fiV{r) + tlj{r} = (4.1) 

where fi is the reduced mass of the particles, V{r) is the potential between the two 
particles and k is the relative momentum between the particles in the absence of a 
potential, given by k"^ = 2fiE, where E is the kinetic energy of the system. Recall 
that we are using the particle physics convention that c = h = 1. 

If we now select a solution with a particular orbital angular momentum and set 
ipim{r) = Ri{r) Yimi^r)^ then Ri{r) satisfies the radial equation 



1 d^ 
r dr^ 



2fiV{r) + fc^ 



Ri(r) = 0. 



(4.2) 



If we further set Ri{r) = yi{r)/r, then yi{r) satisfies the equivalent radial equation 



/(/ + 1) 



dr^ 



2/iV(r) + k'^ 



yi[r) 



(4.3) 



The solution of Eq. |4.1| for V{r) = is just a plane wave, which can be broken 
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up into its partial waves, 

-ik ■ f °° 

cP^ir) ^ e = 4vr^ 5: t^^kr) Y*^^{n,) Y,^,m, (4.4) 

where ji{kr) is a spherical Bessel function and k gives the direction of the incident 
particle. The equivalent solution of Eq. [4.2| for V{r) = (picking the only particular 
solution that is regular at the origin) is 

Ri{r) = ii{kr) ^T^^ — gin {kr — \Iti). (4.5) 

For the case with a potential, the stationary scattering waves ■^/'^ (r) are solutions 



of Eq. [4.1| . Their asymptotic forms consist of the incoming plane wave, and either an 
outgoing (+) or incoming (— ) spherical wave. We will be interested in i'^{r) because 
it will be describing the final state, and hence must be propagated backwards in time 
as an incoming spherical wave to be projected onto the initial state plane waves at an 
equal time. In nonrelativistic collision theory it is defined to have the following 
asymptotic form 

^ —ikr 
^,-(r1™e^^-%/r(a,)^ (4.6) 

where f^{Qr) is called the scattering amplitude and describes the effect of the scat- 
tering. 

By symmetry, ip^{r) and f^{Qr) will be symmetric about the incident axis. It 
is then reasonable to choose to set the incident axis along z (i.e. k = kz), so that 
neither function is dependent on 0^- We can then make the following expansions in 
terms of Legendre polynomials 

1=0 ^ 
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oo 

fui(^r) = Y.frPlicOs9r). (4.7) 
/=0 

where yi{r) is a regular solution of the equivalent radial equation (Eq. ^^ ), and has 
the following (un-normalized) asymptotic form^ 

yi{r)'^ ai sin {kr - ^In + Si). (4.8) 

We can then match the two asymptotic forms of ilj^{r,6r) to obtain ai and ff. 
Then applying the spherical harmonic addition theorem to the Legendre polynomial 
gives for the stationary scattering wave (with incoming spherical wave) solution of 
Eq. 13 

oo I ■ r 

V',-(r1=47r^ Y: ^'e ' ui{k,r)Y;^^{n,)Yi^^inr) (4.9) 

1=0 mi=—l 

where ui{k,r) is the regular solution of Eq. |4.2| with asymptotic form 

ui{k,r)"^ ^sm{kr-U7T + Si). (4.10) 
kr 

4.3 The Fermi Approximation 

The simplest way to perform a correction due to FSI's is to use the Fermi approxi- 
mation We assume that the FSI interaction is very short range, and hence that 
only the value at the origin of the relative wavefunction of the final state particles is 
significant. Since the wavefunction enters into the amplitude, the S-wave cross-section 
would be enhanced by a factor of |uo(^, 0)/jo(0)P = \uo{k, 0)p, and the cross-sections 
of other partial waves would be unaffected. 

^We make the usual assumptions about the potentiah that it is analytic in the vicinity of the 
origin, that lim r'^V{r) = and that lim rV{r) — 0. 

r — ^0 r — *oo 
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4.4 Applying FSI's to the Flux-tube Breaking 
Model of Meson Decay 



decay. Although we do not apply their methodology in this work, we give the details 
here for completeness and because: it connects the two topics of this thesis (meson 
decay and FSI's), we are aware of no detaz/ed published discussion of the method, and 
our derivation includes factors neglected by Geiger and Swanson. 

Recognizing that the last exponential function of Eq. p.23| is the wavefunction 
between the final state particles in the absence of an FSl potential (Eq. [4.4] ), Geiger 
and Swanson replace it with the wavefunction in the presence of such a potential 
(Eq. [4.9|) . By doing so they are using a final state in which the FSI effects are 
included, instead of simple plane waves in which the effects are ignored. This should 
give a final result that includes the FSI effects. 

Let us first examine the partial wave amplitudes from the fiux-tube breaking 
model in the absence of FSI's. In Section |2.4| , we used the Jacob- Wick formula to 
convert to partial-wave amplitudes (see Appendix ^ - here we use the recoupling 
calculation because we have a need for the / dQ Y^^j (fi) integral. This gives us, from 



Geiger and Swanson have applied FSI's to the fiux-tube breaking model of meson 



Eqs. ^ and using Eqs. and [E^ 




(P) = 70 ^SEaEbEc E {LMlSMs\JaMj^) 



Ml A Ms A '^Lg MSb Mjb . 
Ml ,Msr MlMs,'^. 



X {JsMjg JcMjJSMs) {LaMl^SaMsMMj^+Mj^)) 
x{LbMlbSbMsJJbMjJ {LcMLaScMsJJcMj^.) (lml-m|00) 



(4.11) 



Chapter 4. Effects of Final State Interactions 



65 



where 



I^^\P, mi, 1712,1712,) 



(27r 



df dw iJn^LsMri-w - r) i^l^L^uAw ' r) 



2 mm 



iP • (m+r + m_w) 

X c 



[rA) 



(4.12) 



Note that in evaluating M^^{P) we can pick any value of Mj^; ahernatively, we could 
sum over Mj^ and divide by {2 J a + 1), on the right side. 
Next, set 



Pe 



iP ■ (m+r + iri-w) 



i ^ iP ■ {m+rBc + rri-w) 



1714 



rBC=r 



then use an alternate form of Eq. [4.4| to set 



(4.13) 



1=0 mi=—l 



(4.14) 



where s = mj^fBc + ^~w, and finally, apply the orthogonality of spherical harmonics 
to obtain 



/^^^(P, mi,m2,m3) 



(27r 



m. 



w] 



xe 



2""-^mm 



- rA=-2r,rsc=*'/' 

(4.15) 



where, in the absence of FSI's, 



(4.16) 
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Eq. |2.26| can be used to break Ylm{^s) up into spherical harmonics with argu- 



ments ^i'ga and of course, s is given by 



s = + w?_w'^ + 2m j^m-TBcw cos 9 (4-17) 

where 9 is the angle between fsc and w. 

The above is just another way of expressing the flux-tube breaking model of 



meson decay of Section p.3| . However, we can include the effects of FSI's by replacing 
the plane wave ^^P-i'^+'^Bc+m-w) ^\^}^ V'^(^) (Eq. |4.9| ), instead of just expanding it 
with Eq. We then get, instead of Eq. [4.16| ,^ 



^{BOLM^ifBc.w) = 47rr^e' ^ ml(P,s) FlmJ^^s). (4.18) 

In the absence of an FSI potential, 5l =^ 0, ul{P,s) =^ jiiPs), and we recover 
Eq. |16. 



4.5 Applying FSI's to QED Amplitudes 

In order to develop the expression for applying FSI corrections to QED amplitudes, 
we first find the expression for correcting amplitudes in nonrelativistic collision theory 



34| , in the Born approximation. Once the expression for the correction is derived, we 
use the cross-section to relate the collision theory amplitude to the QED amplitude. 
Applying the resulting relation gives us an expression for correcting QED amplitudes 
for FSI effects. We specialize to an interaction with two particles in the final state. 
In collision theory, consider the Hamiltonian H = Hq + V where Hq is the kinetic 



■^Note that Geiger and Swanson associate a different constant factor with their definition 
of iP{bc)lMl{^bC:W)! they neglect the phase and they do not exphcitly mention the factor of 

iSl 

e 
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energy term, and V is the potential. V — W + U is divided into the potential 
representing the basic interaction, W, and the FSI potential U. The plane waves 
0g(r) are the solutions of the Schrodinger equation in the absence of a potential, 

{Ho - E) <l>^,{r) = 0, (4.19) 

the stationary scattering waves X^(0 are the solutions with the full potential, 

{Ho + W + U-E)x^{f^^O, (4.20) 

and the stationary scattering waves i^^if) are the solutions with the FSI potential 
only, 

{Ho + U -E)ij^{r) = Q. (4.21) 
The X^(0 solutions of the equivalent integral equation. 

The regular transition amplitude is defined in terms of the transition operator be- 
tween plane waves, or in terms of the potential between a plane wave and a stationary 
scattering wave: 

(/|TK) ^ (0,-jT|0,^; = {<t>^,^\V\xl) = {xl\V\<t>^,). (4.23) 

In the Born approximation, wc assume that V is sufficiently small that X^if) — 
and we can replace xi(r) by 0^(r) in the expression for (/|T|i): 

(/|T|z) (/|r|i)^ORN ^ ). (4.24) 
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We will use the Bom approximation from this point on, but will drop the ^^rn 
notation for convenience. 

Now consider the problem of FSI's. We have a very good theory (QED) that will 
calculate the transition amplitude due to just the basic interaction between plane 
waves, {f\T\i) = (0r but we need to find a way to include the FSI effects 

(the potential U). To do this we use a final state that includes the FSI effects, V'^(r), 
so only the potential W is left between the initial and final states: 

{fmr'={^^^\W\<f>^,). (4.25) 

Next, we make use of the plane wave orthogonality and closure conditions 



1 



^g) = {2'Kf5\k-k'), (4.26) 



d^A: = / (4.27) 



to obtain 



^ VfcV^r(^)(/|TK), (4.28) 



(2vr) 

where il)^*{k) is the complex conjugate of the momentum-space solution of the Schro- 
dinger equation (Eq. |4.1| ) with the FSI potential U . 

To be able to apply this to QED amplitudes, we need to know how the non- 
relativistic scattering amplitudes relate to QED amplitudes. A comparison of the 
expressions for the total cross-section reveals that |A^p oc s|(/|T|i)p where s is the 
Mandelstam variable.^ Fourier-transforming Eq. [4.9| to momentum space, we find 



*In the CM frame, the total energy is given by y^. 
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that the QED amphtude corrected for FSI effects is 

^ i=omi=-i '' ys{k) 

X drr^ji(kr)ui(kf,r) (4.29) 

where s and A4 have been written as functions of the momenta in order to show 
whether they are calculated with the variable of integration, k, or the relative mo- 
mentum between the final state particles in the CM frame, kf. There is ambiguity 
as to how to relate k and s{k) - nonrelativistically or relativistically. We follow the 
usual practice of keeping as much as possible relativistic, and so use for equal-mass 
final state particles: s{k) = Ak"^ + 4m^. 

It is worthwhile to expand each M. into partial waves, 

oo L 

M{k)^Y. E fLMAs{k))Yr.MA^k) (4.30) 

L=0 Ml=-L 

in order to get the FSI correction for a single partial wave 

flUsikf)) = l^)e^' Hdk Hdrr^e l^^^^^^ j,{kr) u,{k,^r). (4.31) 

TT JQ JO \/s{k) 

In the absence of an FSI potential, Sl =^ 0, UL{kf,r) =^ jL{kfr), the orthogonality 
condition of the spherical Bessel functions gives us a delta function, and we obtain 
fLMr{s{kf)) = fLuMkf)) as we should. 



Chapter 5 

Final State Interactions: An 
Application 



5.1 The Interaction 77 tttt 

In 1986, a group at Orsay (using DM1 on DCI) measured the cross-section of 77 
TT+TT^ near the vr vr threshold, and found a cross-section about twice that expected 
from scalar QED^ Although statistically not very significant (a 2a effect), the 
result aroused interest because it differed from theoretical expectations. By Low's low 



energy theorem the interaction of low energy photons depends only on the static 
properties of the target. Because at low energies the long wavelength photons cannot 
make out the individual constituents of the pions, the pion-photon coupling depends 
only on the charge of the pion, and the pion can be treated as a point particle in 
scalar QED. The strong final state interactions (FSI's) between the outgoing pions 
were not expected to be as significant as the electroweak effects. 



^Scalar QED refers to a subset of QED that deals with the interactions of just photons and scalar 
(spin 0) particles. Because the pions really consist of two spin-^ quarks, scalar QED is used here as 
an effective theory that is valid at low energies because the long wavelength photons cannot make 
out the individual quarks. 
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Figure 5.1: Graph of the 77 tt+tt^ total cross-section vs. ^/s, for the tree-level 
scalar QED and BDI predictions (reproduced by us), and recent data. The curves have 
been corrected to have a limited polar acceptance to match the data: | cos6\ < 0.6. 
The data is from Behrend et al. |^8|, Boyer et al. and Aihara et al. Some 
of the data shown in this chapter was obtained by us from References |^ and 
The horizontal error bars on the data show the bin sizes; for the vertical error bars 
all of the given errors were added together in quadrature. 

In 1987, Barnes, Dooley and Isgur (BDI) ||4^ estimated the effects of these FSI's 
using the Fermi approximation (see Section [4. 3] ) with an effective potential extracted 
from a quark model, and found a similar enhancement. Unfortunately, the experi- 
mental enhancement has since disappeared; the early experiments were apparently 
plagued by problems with subtracting the lepton background due to low statistics 

and the later experiments show no such enhancement. Figure |5.1| shows the two 



curves and the latest experimental data. Note that the data fits the tree-level scalar 
QED curve well for the whole range shown (up to a CM energy of 900 MeV), even 
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though we would not necessarily expect it to be accurate at the higher energies. 

Calculations of the 77 tttt cross-sections have been done using other methods. 
One-loop chiral perturbation theory has produced mixed results: the 77 tt^tt" 
cross-section curve shows a small enhancement over the tree-level scalar QED 
curve near threshold that is still compatible with the latest experimental data, while 
the 77 TT^TT^ cross-section curve ^ has a shape that is rather different from the 
data. However, a two-loop chiral perturbation theory calculation of the 77 ti^tt^ 



cross-section agrees with the data. The difference is apparently mainly due to 
TTTT rescattering and renormalization of the pion decay constant. 

Calculations of the 77 —>■ tttt cross-sections have also been done with dispersion 
relations, typically using phase shift data from tt-tt scattering, and constraints from 
unitarity, analyticity and crossing. The 77 — tt+tt^ cross-section results are again 
quite similar to both the tree-level scalar QED curve and the experimental data near 
threshold The 77 tt^tt^ cross-section results also do well [0, showing 

good agreement with the data. It is believed that higher order exchange effects, 
whose absence is problematic for the one-loop chiral perturbation theory calculation. 



are included in the dispersive calculation through the unitarity requirement [ ^ . 

On the face of it, the failure of the quark model to do as well as these other two 
approaches represents quite a blow to its reputation. Being able to predict the effects 
of FSI's is important if we are to understand QCD. In addition, if we understood the 
situation in 77 two pseudoscalar mesons, we could then approach the situation 
in 77 two vector mesons with more confidence; the structures seen in the cross- 
sections of these processes are not well understood and are a long-standing puzzle. 

However, there is room for improvement in the BDI calculation. We will redo the 



calculation using the full FSI apparatus of Section ^75| , and newer, better potentials. 
In addition, we will calculate the cross-section for the process 77 — tt^tt^. Since 
neutral pions do not couple to long wavelength photons (there is no such process in 
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(a) 



(b) 



(c) 



Figure 5.2: The tree-level Feynman diagrams for the scalar QED interaction 
77 — > '7r"'"7r~. The labels of the particle momenta (pi) and photon polarization vectors 
(ej) are shown. 

scalar QED), 77 — > tt^tt^ is expected to be suppressed relative to 77 tt+tt^ making 
the effects of the FSI's relatively much more significant. In our calculation, it is the 
FSI potentials that mix the tt+tt^ final state into tt^tt^. 

TT+TT^ in tree-level scalar 



In Section |5.2| we calculate the cross-section for 77 
QED. In Section ^]3| we discuss the effective potentials and the application of the full 
FSI apparatus (see Section ^TS]) to this problem. In Section we give the results of 
our calculations. 



5.2 77 TT+TT in Tree-level Scalar QED 



The three tree-level scalar QED Feynman diagrams^ contributing to the interaction 
77 —>■ TT+TT^ are shown in Figure |5.2| . 

The corresponding Feynman amplitude (using the momenta and photon polar- 
ization labels shown in Figure |5.2|) is 



M = 47ra 



ei ■ (2p3 - pi) €2 ■ (2p4 - P2) ei ■ (2^4 - Pi) £2 ■ {2p3 - P2) , ^ 

H ; 7. r /Ci ■ 62 



{P2 - PiY - rn^ 



{pi — PaY — 'm^ 



(5.1) 



See Appendix A. 2 for a brief explanation of Feynman diagrams. 
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Figure 5.3: The geometry of the 77 — tt+tt interaction in the CM frame. 

where a is the fine structure constant and m is the pion mass. Averaging the photon 
polarizations and integrating over phase space gives us the total cross-section 



a 



TXOL 

4m^ 



— X — X (2 — x) In 



'1 + 



X 



\-^f\ 



X 



(5.2) 



where x = and s is the Mandelstam variable. 

s 

We would now like to work out the cross-sections corresponding to different partial 
waves in the final state. In the CM frame, we take pi to point in the positive z 
direction, and the direction of ps to be described by the spherical polar coordinates 
6 and 0, as shown in Figure |5.3|. For the photons we take circular polarizations: 



ei(+)=e2(-) = -^(0,-1,-^,0), 
6i(-) = e2(+) = -^(0,1,-2,0). 



(5.3) 



We calculate the helicity amplitudes 

Sua X 



1 1 



1 — (1 — x) cos^ 6 

±2u 



, , Sua (1 — x) sin^ 9 e~ , , 
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and then expand them in a spherical harmonic basis to obtain the partial waves (e.g. 
for the ++ case) , 



M 



L=0 Ml=-L 



(5.5) 



We find that fj^^ = fj^j^ is only non-zero for = and that / 



LMr 



fl-Mr 



only non-zero for = 2.^ The cross-section is then written in terms of the partial 
waves as 



a — 



\/l — X 
1287r2s 



(/oVr+ E [(/.Vr + (/.V)^ 



L>2,even 

The needed partial waves (see below) are given by 



(5.6) 



++ 

00 



f++ 

J 2 



120 



22 



f++ 
JAQ 



42 



/■++ 
J 60 



3 

87r2Q; 



X 



In 



X 



1 + vr^ \ 



4-\/57r2Q; ^ 



A 15 3 
4W —'K'^a 



_ 3 X 



1 — X 
10 



2 + x 1 + 
— D H , m 



\/l — X ^ 



X 



X 



1 — X 



3V107r2a 



„ /— 3 X 



3 1 
110 

~3~ ~ l-x 
-54 76 



{l~x) 



In 



1 + 



a; 



1 - Vl^. 



70 3a;2 + 24a; + 8 , / 1 + ^/l^' 
+ 5 In ' 



14 x'^{<6 + x) 



In 



1 + 



X 



l-x 



3(1 -x) (l-x)2 (l-a;)l VI "V^ 
231 7(4 -15a;) 21(1 + 5a; + 5a;2) 



X , 



40 



+ 



8(1 - a;) 



^{l-xf 



16 + 120a; + 90a;2 + 5a;3 A + ^/T^' 
+ : In' 



(5.7) 



16(1 -a;)t 

The total cross-section and the cross-sections of the L — 0,2,4 partial waves are 



■^The orbital angular momentum L between the pions is even because of charge conjugation 
invariance. 
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Figure 5.4: Graph of the 77 —>■ tt^tt^ tree-level scalar QED cross-sections vs. ^/s. 
The total cross-section is shown, as well as the cross-sections of the L = 0, 2, 4 partial 
waves, and the cross-section corresponding to the sum of the L = 0, 2, 4 partial waves. 
The curves are for full polar acceptance: | cos^| < 1, which is why they differ from 
those shown in Figure 



shown in Figure 5.4. 



5.3 Applying the FSI Correction 



5.3.1 The Effective Potentials We Use 



The 71-71 potentials that we use were calculated by Swanson et al. |52[ using 

techniques developed by Barnes and Swanson . The potentials come from equating 
the T-matrix giving the leading order interactions in a perturbation theory for a 
nonrelativistic quark model, to the T-matrix for the interactions of point-like mesons. 
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c 



D 



a) 



b) 



Figure 5.5: Example diagrams of the two classes of interactions included in the quark 
model used to find the effective potentials: a) t-channel gluon exchange, b) s-channel 
gluon exchange. 

written in terms of the effective potential they were trying to find. 

The potentials are different for the different isospin states of the two pions. Be- 
cause of charge conjugation invariance, the orbital angular momentum L between 
the two pions must be even. Because of Bose-Einstein symmetry (which holds for 
members of the same isospin multiplet in the context of the strong interaction ||54||), 
the sum L + I must be even, and hence the total isospin I must be even as well. It 
then follows that the isospin states (written as I//2)) of the pions are 



|7r'''7r ) 



--L|00) + /||20>. 



(5. 



There are two classes of interaction included in the quark model used to find the 
effective potentials. The first, referred to as t-channel gluon exchange by Swanson 
et al, involves both one-gluon-exchange and a linear confining potential, followed by 
quark rearrangement back into colour singlet states. This can occur for both / = 
and 1 = 2 states. The second, referred to as s-channel gluon exchange, involves the 
annihilation of two of the quarks, producing a meson hybrid in the intermediate state. 
Because of the intermediate state, this interaction can only occur for 1 = states. 



Example diagrams of the two classes of interaction are shown in Figure 5.5. 
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We only present results for FSI corrections to the L = partial wave of 77 vrvr 
because the L = 2 potentials available from the same sources are smaller than those 
for L = 0, and we found that their effect on the cross-section was negligible. 

Although relativistic phase space is included in the quark model used to find the 
effective potentials, the potentials themselves are essentially nonrelativistic. It is our 
hope that they will contain enough physics to correctly calculate the distortion of 
the final state wavefunction leading to FSI effects. However, it is known []55| , |53| , pO| 
that the t-channel gluon exchange potential is unable to accurately predict the phase 
shifts resulting from 1 = 2 tt-tt scattering, apparently because it lacks relativistic 
phase space. It is expected that the other potentials suffer from the same problem. 

Weinstein and Isgur have argued that this is related to the inability of the 



nonrelativistic quark model to accurately predict the charge radii of the ground state 
mesons and baryons, and so apply a correction factor n = 2.0 to the ranges of the tt-tt 
potentials, and an overall normalization factor rj = 0.6 which is fitted to experiment. 
To be specific, a potential originally parametrized by 



V{r) = Voe H'-o; (5.9) 

is written as ^ 

V{r) =T]Voe'^(^) (5.10) 

with the corrections. However, it is not clear whether a factor that corrects the phase 
shift for the nonrelativistic phase space will necessarily give the correct dynamics for 
the distortion of the wavefunction. A better approach might be to use the potentials 
without T] and k to obtain the wavef unctions, but to use the phase shifts calculated 
directly in the quark model^ to predict the FSI effects (in Eq. [4.31| for example). 



Swanson et al. M, BO, 51 1 found that the quark model, with its relativistic phase 



We refer to the specific quark model used by Swanson et al. to extract the effective potentials. 
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Figure 5.6: Graphs of the tt-tt potentials used in this work vs. r. Note the different 
scales on the two r axes. T-channel gluon exchange is present in all potentials. 
For / = 0, the potentials are plotted both with and without the contribution from 
s-channel gluon exchange. In a), the potentials without the factors 77 and k are 
plotted. In b), the potentials with rj — 0.6 and k, — 2.0 are plotted, except for one 
curve where the contribution from t-channel gluon exchange contains 77 and k, but 
that from s-channel gluon exchange does not. 

space, does accurately reproduce the tt-tt scattering phase shifts away from threshold 
(although resonance effects must be added for / = - see below). 

In addition, since 77 and k have previously been applied only to potentials rep- 
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71-71 Potentials [Vq in 


GeV, 


ro in GeV-^) 


Pion State 


t-channel gluon exchange 


s-channel 






hyperfine 


confinement 


gluon exchange 


L 


I 


Vo 


















-0.392 


1.36 


-0.024 


2.29 


-1.75 


1.48 





2 


0.786 


1.36 


0.047 


2.29 






2 





-0.044 


1.40 


0.0175 


1.49 






2 


2 


0.088 


1.40 


-0.035 


1.49 







Table 5.1: The parameters of the potentials used in this work [S^, ^ The 
t-channel gluon exchange potentials have two contributions, due to colour-hyperfine 
and confinement terms, which must be summed. The potentials are parametrized as 
in Eqs. ^]9| and 5.10| , depending on whether rj and k are to be included or not. 



resenting t-channel gluon exchange, it is not entirely clear if they should be applied 
to the potential representing s-channel gluon exchange. We will examine these three 
possibilities below. The various potentials, with and without the t] and k, are shown 
in Figure |5.6| . The parameters Vq and ro of the potentials considered (including the 
L = 2 potentials, for completeness) are given in Table |5T 



5.3.2 The FSI Apparatus Applied to the Pions 

Calculating the effects of FSI's on the 77 tt+tt^ cross-sections is complicated by 
the fact that there are two isospin states involved. Because the / = and 1 = 2 
states feel different potentials, they will be affected differently by the FSI's. This will 
unbalance the isospin combination that is tt+tt", leading to some Tr^7c^ production. 
BDI used a single tt+tt^ potential, and hence neglected conversion to tt^tt^. 

We have two choices as to how to proceed. We could transform the potentials from 
the / = and 1 = 2 basis to the vr+vr^ and 7T^7t^ basis (leading to off-diagonal potential 
terms that would mix the two states), and then solve a coupled-channel Schrodinger 
equation to find the final state wavefunctions needed for the FSI calculations (see for 



example Reference [0). Alternatively, we could take advantage of the facts that the 
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potential matrix is diagonal in the isospin basis (since isospin is conserved in strong 
interactions), and that we are assuming the pion states to be degenerate in mass (at 
137 MeV), and just solve uncoupled Schrodinger equations in the isospin basis. We 
opt for this second choice. 

The cross-section for 77 — > tttt, corrected for FSI effects, is given by ^ {of. Eq. 5.6 ) 



a: 



FSI 



1287r2s 



fFSI 

J 00 TTTT 



+ 



00 



L>2,even 



fFSI 

J L0-K1T 



+ 



fFSI 

J L2 -KIT 



(5.11) 



Note that we have not yet specified the charges of the pions - S = 1/(1 + ^tttt) is a 
statistical factor which is needed if the pions are identical particles (i.e. ti'^tt^). 

The calculation of flfj^ proceeds similarly to that in Section ^]5| - the first 
divergence is in Eq. ^.28| , which is replaced by 



FBI 



(27r)3 

^/d3M^-|0f^-) 



) 



(5.12) 



where we have inserted plane waves for both tt+tt" and -n^ix^ - the tt^tt" term is then 
dropped because vr^vr^ isn't produced from 77 in tree-level scalar QED, which is what 
the potential W represents. 

We can proceed as before, and relate the final bra-ket to the scalar QED ampli- 
tude, and the left-hand side to the corrected scalar QED amplitude. The first bra-ket 
can be written in the isospin basis using Eq. |5]3 as one of 



) 



'2 ^^-0| ^ 1 ^^-2| 

3 ' 73 



'^Note that we have dropped the and + hchcity notation on the /'s - it can be identified 
from the value of Ml anyway (A/l = ^++, Ml = '2. ^+"). 
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(27r 
(2vr: 



3 ^fc/ ^ ^ 3 ^' 



(k) 



(5.13) 



depending on which pion state is produced in the end. Here i/j^; *{k) is the complex 



conjugate of the momentum-space solution of the Schrodinger equation (Eq. |4.1| ) 
with the potential for isospin /. If Ui{kf,r) is the solution of the equivalent radial 
equation (Eq. [4.2|) with orbital angular momentum L, then Eq. [4.31| becomes for 
TT+TT^ production, 



'} I roo roo 



sik) 



jiikr) 



X 



2 i6 
- e " -u 



L „,o 



1 iS 



ri + - e ^ 



{kf,r) 



where 5£ is the phase shift of u^j^ikf^r). 
In order to keep things real, we define 

/I U I' u 



sik) 



which gives 



^FSI 

J LMl tt+tt- 



2 _ 4 
~ 9 



{alM^f + I {alMr)^ + H ^iA/L ^LAf^ COS (5^ - 6 



(5.14) 



jLikr)u[{kf,r), (5.15) 



(5.16) 



Similarly, for vr vr production, we get 



fFSI 

JLMLTr°TT° 



2 2/q \2 2/2 \2 4q 2 /e-n e-2\ 

(^lmJ + q (^lmJ - q ^lm^ 9lMi^ cos (5^ - 5^). 



9 



(5.17) 



Note that if there is no FSI potential, (5£ =^ 0, u^i^{kf,r) => iiikfr), and the 
orthogonality condition of the spherical Bessel functions gives us a delta function, so 
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9lml =^ fLAiL and we obtain = JlMl and flfj^^o^o = 0, as expected. 

In this work, we only correct the L = partial wave, because higher partial waves 
have negligible FSI effects. Thus in Eq. |5.11| only the /qqI^, term survives for 7r°7r°, 
and for tt^tt^ it is the only term that differs from the uncorrected value. 



5.3.3 The Problem of Limited Polar Acceptance 

Expanding the scalar QED amplitudes in a spherical harmonic basis (Eq. ^.5|) to 
obtain cross-sections for particular partial waves is essential to our calculation, but it 
requires that we be able to integrate over the entire An angular range when calculating 
the total cross-section. Unfortunately, experiments don't measure data over the whole 
polar angle (| cos^^| < 0.6 — 0.8 is typical), so we cannot easily compare our results 
with experiment. 

Integrating to get the total uncorrected cross-section for a limited polar accep- 
tance (— cos6'acc < cos 6* < cos6'acc) is uot a problem: Eq. |5]^ is replaced by 



2 cos 6, 



Am? 

— x^(2 — x) In 



acc Xy\ X 

1 + 



X 



1 - (1 -X)C0S2^« 



X COS Of. 



1 — ^/l — X COS 9a 



(5.18) 



However, the expression in terms of partial waves is affected more drastically. 
Eq. |5.6| becomes (where once again we drop the explicit helicity labels) 



0"a 



2tt rcos 
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(5.19) 



For the limited polar acceptance, the spherical harmonics are not orthogonal to each 
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other, so each magnitude contains an infinite number of terms. As well, the contri- 
butions of the partial waves to the total cross-section can no longer be separated. 

Because the L = partial wave is expected to be the only one affected significantly 
by FSI's, we assume that all of the observed 77 —* tt^tt^ events are in an L = state. 
Since that distribution is spherically symmetric, we can correct the data to full polar 
acceptance by simply dividing it by the value of cos6'acc appropriate for the particular 
experiment. 

For 77 TT+TT^ the total cross-section observed is a mixture of all partial waves; 
again, we assume that only the L = wave is affected by FSI's. Define 



h] 



2n 



cos 0a 



COS 9s 



d(cos^) Y*, 



(5.20) 



and note that hf^i^ = h^i^ and /i^"^^ = ^lih- ""^^ then write the total cross-section 
with only the L = partial wave corrected for FSI's as the total uncorrected cross- 
section, plus an infinite number of correction terms, each of which is the difference 
between the corrected and uncorrected values of a term in Eq. |5.19| . Only those terms 
involving /oo, which goes to /^^^+^_ in the correction, need correcting. Eq. ^.11 



IS 



replaced by 



a. 



FSI 



'^acc ~l~ 



"■oo 



fml+TT- ~ (/oo)' 



°° * 



L>2,even 



4 2 X 24 
hm\g{9m) + g (^?oo) + g ^?oo ^?oo cos (5° - ^o') 



(/oo)' 



°° (2 1 

+ E 2 ho \ ^ 9oo cos 5° + - ^^0 cos Si - /oo 

L>2,even ^ 



(5.21) 



We find that we get excellent convergence over the energy range we are interested in 
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when we truncate the series after L = 6. Fortunately, the three experiments whose 
data we use for comparisons all have the same limited polar acceptance (| cos < 0.6), 
so we only have to correct our results for one value of cos6'acc- 

5.3.4 Numerical Details 

We solved the radial Schrodinger equation (Eq. [4.3|) with the appropriate potentials 
using the Bulirsch-Stoer method for differential equations (with Stoermer's rule for 
2nd Qj-(^Qj- conservative equations) |^ . The solutions were started off from the origin 
with the forms of r times the spherical Bessel functions of the same L, since the effects 
of the potentials are insignificant at the origin. The amplitudes and phase shifts of 
the solutions were extracted by fitting the last oscillatory cycle to a general sinusoidal 
function. We carried out the integrals over r and k using an extended Simpson's rule. 



5.4 Results 

We have calculated the effects of the FSI's on 77 tttt for three primary cases of the 



FSI potential (see Section |5.3.1| for a full discussion). The three cases are: 



Case A. In this case we do not use 77 and n in the potentials. We know that the 
potentials will then fail to accurately predict the phase shifts. In order to avoid 
problems from this, we use the phase shifts directly calculated in the quark 
model in our expressions for FSI corrections. For more information, see the 
detailed discussion on rj and n on pages 78 through 80 . 

Case B. In this case we use rj = 0.6 and k = 2.0 in all of the potentials, and we use 
the resulting phase shifts in our expressions for FSI corrections. 

Case Bl. This case is the same as case B, except that rj and k, are used in the 
t-channel gluon exchange contribution to the potentials, but not in the contri- 
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bution from s-channel gluon exchange. 

In addition, we have three secondary cases, A', B' and Bl'. Each of these is similar 
to the corresponding unprimed case, but does not use the same phase shifts in our 
expressions for FSI corrections. Instead, they use a simple (straight-line) expression® 
for the phase shifts that attempts to describe the experimental data. We do this in 
order to get an idea of how much our results are affected by inaccurate phase shifts. If 
our results are poor because we have ignored a particular contribution to the potential 
(such as resonance effects), then using accurate phase shifts in the expressions for FSI 
corrections would partly correct for this. Of course, we would still need to determine 
the effect of the neglected contribution on the distortion of the wavefunction, and 
hence on the FSI effects. 

Let us first examine the phase shift results. The 1 = 2 phase shifts are shown in 



Figure |5.7| . For the three cases A, B and Bl, the curves shown are the phase shifts 
that come from solving the Schrodinger equation with the potentials. The phase shift 
for case A is poor, as expected; in the expressions for FSI corrections we will replace 
it with the curve calculated directly in the quark model. The phase shifts for cases B 
and Bl are identical because they only differ in whether or not rj and k are included in 
the contribution to the potential from s-channel gluon exchange, which is not present 
for 1 = 2. For the range of data shown, it is not clear which is better: the B and Bl 
phase shifts from the Schrodinger equation, or that from the quark model. However, 
data at somewhat higher energies favours the quark model curve. The remaining 
curve shown is just a straight-line which we used to describe the data for the primed 
cases. 



^The straight-line expressions we use to describe the phase shifts are: 
1=0: 5g = 0.0027 rad/MeV x (^i- 250 MeV) 
1=2: Si = -0.00062 rad/MeV x (^i- 250 MeV). 



Chapter 5. Final State Interactions: An Application 



87 



Case A (calculated) 
Cases B, B1 (calculated) 
Quark model prediction 
Straight— line approximation to data 
Hoogland 77 (method A) 
Hoogland 77 (method B) 
Prukop 74 



^-2.5 - 

-5 ^ 

-7.5 r 

-10 ^ 

-12.5 r 

-15 ^ 

-17.5 r 

-20 - 

□ 

-22.5 ' ° 

A 

_25 t \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 

200 300 400 500 600 700 800 900 

Vs (MeV) 

Figure 5.7: Graphs of the 1 = 2 tt-tt scattering phase shift vs. ^/s, for our predictions, 
the quark model prediction, and experimental data. The data is from Hoogland et 
al. |^8| - we show the results of both of their methods; and Prukop et al. []59 
show the results of their first fit. For additional comments see Figure |5.il| . 



- we 



The / = phase shifts are shown in Figure |5.4 Note the rather unfortunate 
behaviour of the phase shift for case B, which suggests that t] and n should not be 
used in the s-channel gluon exchange contribution to the potential. Since there is no 
a priori good reason for not using them in this part of the potential, we must question 
their use at all. However, we will continue to show the results arising from cases B 
and Bl. Again, the phase shift for case A is poor, and for FSI corrections we will 
replace it with the curve arising directly from the quark model. Note that this time 
the quark model curve is low as well; the reason is as follows. When the quark model 



predictions for the / = phase shifts were calculated three contributions were 
considered: t-channel gluon exchange, s-channel gluon exchange, and contributions 
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Figure 5.8: Graphs of the / = n-n scattering phase shift vs. ^/s, for our predictions, 
the quark model prediction, and experimental data. The data is from Mukhin et al. 
|60| ; Rosselet et al. |^ - the horizontal bars only approximate their bins, and their 



data is actually for 5^ — 5\ - we have used the 5^ data extracted from it by Li et al. 
51| ; Estabrooks and Martin [^] - we show the results of both their s- and t-channel 
fits; and Protopopescu et al. - we show the results of their case 1. For additional 
comments see Figure p?T| . 



from interactions with the /o(980) and /o(1300) resonances. These contributions were 
fit to the data, but because we are not including resonance effects in this work, we 
only include the contributions from the s- and t-channel gluon exchange in our quark 
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Figure 5.9: Graph of the 77 —>■ vr+vr" total cross-section vs. y/s, for our predictions, 
the tree-level scalar QED prediction, and experimental data. The curves have been 
corrected to have a limited polar acceptance to match the data: | cos6'| < 0.6. For 



references to the data and additional comments see Figure 5.1 



model curve. When resonance effects are included, the quark model curve fits the 
data rather well. Again, the remaining curve shown is just a straight-line which we 
have used to describe the data for the primed cases. 

Now let us look at the 77 —>■ tt+tt^ cross-section results, shown in Figure [5.9| . 
Cases A, A', Bl and Bl' all fit the data well (as does the tree-level scalar QED 
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Figure 5.10: Graphs of the 77 — > n^n^ cross-section vs. ^/s, for our predictions (cases 
A, A', Bl and Bl') and experimental data. The data is from Bienlein [et al. ) ||64|| , 



Marsiske et al. ||63] and Edwards et al. p6[ . The data has been corrected to full polar 



acceptance: | cos^^| < 1.0. For additional comments see Figure 5.1 



prediction). The effects of using the quark model phase shifts are not significant for 
these cases. Better data would be needed to differentiate between case A, case Bl 
and the tree-level scalar QED prediction. Cases B and B' are disfavoured by the data, 
as might be expected from the poor performance of case B in predicting the phase 
shifts. Because that prediction differed drastically from the data, cases B and B' do 
differ appreciably in their predictions of the cross-section. 

The 77 — > Tx^-K^ cross-section results are shown in Figures |5.10| and ^.11| . Here 
our results do not do as well. This is partly because our 77 — * -k^-k^ cross-section 
is wholly due to FSI's while in the 77 tt^ti^ case the FSI effects are corrections 
to a more significant curve. This also means that the primed and unprimed cases 
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Figure 5.11: Graphs of the 77 ti^ti^ cross-section vs. y/s, for our predictions (cases 
B and B') and experimental data. The upper and lower plots differ only in their 
vertical scales. The data has been corrected to full polar acceptance: | cos6'| < 1.0. 



For references to the data see Figure 5.10. For additional comments see Figure 5.1 



differ substantially. None of our cases are consistent with the data. Cases B and B' 
do particularly badly as expected from our previous results. Of the others, case A' 
comes closest, being roughly a factor of 2 too low. It also shows a local maximum in 
the cross-section at higher energies, in keeping with the data. It is interesting to note 
that using the experimental phase shifts instead of the quark model predictions in 
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our expressions for the FSI corrections greatly improves curve A. Since the difference 
in these phase shifts is due to resonance effects, we tend to beheve that including the 
effects of resonances in the wavefunction distortion as well would improve the curve 
still further. 

We have demonstrated that the quark model can predict the 77 — > tt'^tt' cross- 
section in agreement with experimental data when a sufficiently detailed calculation 
is performed. However, with the present potentials it is unable to accurately predict 
the 77 ^ 7r°7r° cross-section. We find that using Weinstein and Isgur's 77 and k 
factors in the potentials is not helpful when calculating FSI effects; it is better to 
use the expressions for the phase shifts calculated directly in the quark model in the 
expressions for the FSI corrections. We beheve that the next logical step in this 
program of work is to include potentials allowing the I — tttt state to convert to 
the /o(980) and /o(1300) resonances (and back). This would mean solving a coupled- 
channel Schrodinger equation in the / = channel. 



Chapter 6 
Conclusions 



In this thesis we have examined three problems, spht between two general areas of 
the quark model. 

Both of the first two problems involved two models of meson decay: the 3Pq 
model and the fiux-tube breaking model. We have developed general routines for 
these models that can be used for any OZl-allowed strong decay of a meson into 
two other mesons, as long as the radial portion of the meson wavefunctions can be 
expanded in terms of SHO radial wavefunctions, and sufficient computer resources 
are available. The general nature of these routines make them useful tools in the 
study of meson phenomenology. 

While fitting the free parameter of the models and investigating a number of 
choices that must be made in their application, we evaluated the models by com- 
paring their predictions for the decay widths of 28 of the best known meson decays 
to the experimental values. We found that these models are not very accurate - 
the best they can hope for is to predict a decay width to within a factor of 2, and 
even larger deviations are common. Since they are coarse models of a complicated 
theory, this is not surprising. We also found that both the use of SHO wavefunc- 
tions with the effective /?'s of Kokoski and Isgur [Q, and the use of nonrelativistic 
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phase space/normalization give less accurate results than can be obtained with other 
choices. 

The first problem that we investigated was the identity of the /4(2220). Although 
tentatively identified as the '^^4 ss meson by the Particle Data Group, its identity 
has been uncertain since its discovery in 1983. We carried out detailed calculations 
of the decays of the and ^F^ ss mesons in order to see if a meson identification 
for the /4(2220) is tenable. Despite the uncertainties of the decay models, we found 
that the /4(2220), with a measured decay width of approximately 30 MeV, cannot be 
the ss, whose width we expect to be >, 400 MeV. We expect the width of the 
^F4 ss meson to be >, 140 MeV and <, 600 MeV, but the uncertainty of the models 
precludes us from ruling out the identification of the /4(2220) as the ^^4 ss with any 
surety. We do feel it is unlikely however, and propose the following explanation: that 
the broad state seen in hadron beam experiments is the ^^4 ss meson, and the narrow 
state seen in J/ip radiative decay is a glueball. Further experimental data is needed 
to finally identify the /4(2220). 

The second problem that we investigated was the determination of the mixing 
angle between the iri(1270) and ^1(1400) mesons. This was done by comparing pre- 
dictions of the meson decay models for five partial decay widths and two ratios of D 
to S amplitudes with experimental data. We found the mixing angle 6k to be approx- 
imately 44°. If the mixing is due only to the spin-orbit interactions already included 
in the quark model of Godfrey and Isgur, then this implies that the responsible term 
in the Hamiltonian, {Hg^~), has a value far larger than calculated. However, because 
of the delicate cancellation that leads to the small calculated value, we are unable to 
draw any conclusions about whether the mixing is largely due to another mechanism, 
or the quark model Hamiltonian just needs retuning somewhat. 

The third problem that we investigated was the effect of the strong final state 
interactions (FSI's) in the reaction 77 — > tttt near threshold. Using effective potentials 
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extracted from a quark model, we found that our preferred set of potentials (without 
the T) and k correction factors) successfully described the experimental data for the 
77 — > tt'^tt" cross-section. However, the prediction in the absence of FSI's also 
described the data, so better data is need to distinguish between the two. This still 
represents a vindication for the use of effective potentials extracted from a quark 
model, which has been in disrepute because the latest data does not agree with the 
results of a previous, less-comprehensive calculation. We were also able to conclude 
that the use of Weinstein and Isgur's r] and k correction factors for the potentials is 
inappropriate for the calculation of FSI effects. 

Unfortunately, we were unable to reproduce the data for the 77 — > 7r°7r° cross- 
section. This is partly because the 77 — > 7r°7r° cross-section is wholly due to FSI's 
while in the 77 — > tt+tt" case the FSI effects are corrections to a more significant 
quantity. However, we have indications that including the effects of resonance pro- 
duction in the / = channel may bring our results in line with the data. The next 
step in our program would be to solve a coupled-channel Schrodinger equation to 
include these effects. In hght of our poor results for the 77 — > 7r°7r° cross-section, it 
would be premature to investigate the structures seen in the cross-sections of 77 
two vector mesons at this time. 



Appendix A 

Some Tools of Particle Physics 



In this appendix we introduce some of the terms and tools of particle physics, to help 
specialists in other fields in understanding the work of this thesis. For an excellent 
introduction to the field, see Griffiths [E^ . 



A.l Cross-section and Width 

In most particle physics experiments two particles are collided and the outgoing par- 
ticles are collected and identified. Due to the interactions between the particles, they 
may be scattered in different directions, and/or they may be transformed into other 
particles. The probability of a particular reaction occurring is related to the cross- 
section a which is a function of the total energy of the particles and is measured in 
units of barns, or b (1 b = 10~^^ m^). 

A particle will decay eventually if there is an interaction that permits it to be 
transformed into other particles with less total mass. The average lifetime of the 
particle is denoted by r; by Heisenberg's uncertainty principle, this finite lifetime 
leads to an uncertainty in the mass of the particle. In a collision of two particles, 
the production of a short-lived particle in the intermediate state (which then decays 



96 



Appendix A. Some Tools of Particle Physics 



97 



to give the final state particles) shows up in the cross-section plot as a resonance 



described by a Breit-Wigner curve with width F = 1/r (see for example Figure |3.7| ). 

For a particle that can decay into two or more different final states we can define 
its branching ratios, which are the probabilities that the particle will decay into the 
various final states. From these we can define the partial width for a particular decay 
mode as the product of the branching ratio for that mode, and the total width. The 
sum of all the partial widths gives the total width, which is what is used to get the 
particle's average lifetime. We may sometimes refer to a partial width as the width 
for a particular decay mode. 

Similarly for the cross-section for a particular initial state, there are the cross- 
sections going to specific final states, which sum to give the total cross-section. Note 
that in a plot of the cross-section to a specific final state, a resonance in the interme- 
diate state still produces a Breit-Wigner curve with a width equal to the total width 
for that resonance. (Even though we are only looking at a particular decay mode, the 
lifetime of the resonance is still that determined by all of the possible decay modes.) 



A. 2 Feynman Diagrams 

In a perturbation theory such as that normally used to calculate quantities in QED, 
each order of perturbation is represented by Feynman diagrams, such as those shown 
in Figure ^]2| for 77 tt+tt^ in scalar QED. In a Feynman diagram the particles 
are represented by lines, and basic interactions by vertices where a number of lines 
meet. Since for QED each vertex corresponds to another factor of ^/a in the Feynman 
amplitude^ the order of the perturbation is given by the number of vertices in the 
diagram. 

^The Feynman amplitude is similar to the amplitude of quantum mechanics. Its magnitude 
squared is combined with the phase space of the final state particles to obtain the cross-section or 
width. 
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A tree-level diagram is one in which there is only one path of lines connecting 
any two vertices (i.e. there are no loops). For QED at least, if a tree-level diagram 
exists for a particular process then it is the lowest order diagram. 

A. 3 Clebsch-Gordan Coefficients and the Wigner 
nj Symbols 

An angular momentum state, be it spin, orbital angular momentum, or a vector sum of 
some combination of these, is represented (in an SU(2) algebra) by the usual quantum 



numbers \jm) where the magnitude of the total angular momentum is y i(j + 1) and 

the 2;-componcnt is m. To combine two angular momentum states we use the Clebsch- 
Gordan coefficients, defined by 



where {jimij2m2\JM) is the Clebsch-Gordan coefficient and M = mi + m^- 
The Wigner 3j symbol is related to the Clebsch-Gordan coefficient by 






(A.l) 



3=\h-32\ 




(A.2) 



The Wigner 6j symbol 
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is used to combine three angular momentum states, and the Wigner 9j symbol 



jl 32 Jl2 

is j4 Jsi 

Jl3 J2A J 

is used to combine four angular momentum states. We will not define these here, 



though the 9j symbol is defined in Appendix p.3| . For more information, see for 
example Reference . 



A. 4 Conservation Laws and Invariance Principles 

Because of symmetries in the theories of the interactions, the interactions themselves 
obey a number of conservation laws and invariance principles (see for example Refer- 
ence |54|). We will only discuss the ones of interest here. 

There is the usual conservation of energy, momentum, angular momentum and 
charge. 

If one thinks of an antiparticle as being —1 particles, then the numbers of lep- 
tons and quarks are separately conserved as well - this is enforced by the allowed 
interaction vertices (e.g. if an e^e~ pair is created, the change in lepton number is 
(-l) + (+l) = 0). 

Because the strong interaction is independent of the type of quark involved, and 
the u and d quarks have very similar masses, we can think of them as different 
projections of the same state which we represent by the angular momentum-like 
quantity called isospin (where u = |||) and = || —\))- Isospin is conserved in the 
strong interaction to the extent that the masses of the u and d quarks are degenerate. 

The spin-statistics theorem governs how systems of identical particles behave 
when two of the particles are exchanged (boson wavefunctions are unaffected, fermion 
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wavefunctions pick up a factor of —1). It is important that the wavefunctions we 
construct obey these symmetries. In the context of the strong interactions, members 
of the same isospin multiplet also obey the spin-statistics theorem, even though they 
are not strictly identical. 

The parity of a system indicates how the wavefunction behaves under the space 
inversion r — > — r; the wavefunction is either unchanged (T' = +1) or picks up a factor 
of —1 {V — —1). The total parity of a system is conserved in the electromagnetic 
and strong interactions, but not the weak interaction. 

Under the charge conjugation operator, particles are changed to their antiparticles 
(and vice versa). Systems which are eigenstates of C (such as 7r°, or the tt+tt" system) 
may have C = +1 or — 1. Charge conjugation is also conserved in the electromagnetic 
and strong interactions, but not the weak interaction. Charge conjugation can be 
extended by combining it with isospin to get G-parity, which is conserved in the 
strong interaction only. 



A. 5 Relative Coordinates 

The wavefunction of two particles in position-space depends on the coordinates of 
the two particles, fi and r2. When solving the Schrodinger equation it is convenient 
to separate it into equations involving the centre of mass (CM) coordinate of the 
particles, R, and their relative coordinate r: 



R — ri — R + 



mi + 1712 mi + m2 3^ 

^ ^ ^ 5 mi ^ 

r — r\ — r2 r2 — R r 

mi -|- m2 
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The conjugate momenta pi and p2, and P and p, are related by: 

mi 

P = P1+P2 Pi = : P + P 

mi + m2 _ 

m2Pi - mip2 ^ m2 - ^ 

P = \ P2 = ■ P-P 

mi + m2 mi + m2 

This is a nonrelativistic construction: p remains unchanged when both particles are 
given the same velocity boost only for nonrelativistic velocities. 

The CM frame refers to the inertial reference frame in which P = 0; i.e. in which 
the centre of mass of the system is motionless. 



A.6 The OZI Rule 

Okubo, Zweig and lizuka (OZI) noticed that QCD interactions in which all of the 
initial quarks are annihilated, and are not present in the final state, are suppressed. 
This is thought to be because for interactions of this type, asymptotic freedom means 
that the strong coupling constant (which enters the cross-section as at least a^) 
becomes small since all of the energy of the initial state particles must be included in 
the momentum transferred through the gluons (see for example Reference [|6^]). 



A. 7 Meson Mixings 

Particles with the same "good" quantum numbers (i.e. a complete set of commut- 
ing observables) "mix" if there exists an interaction through which one particle can 
be transformed into the other, and vice versa. For example, the strange mesons 
{—us, —ds, —sd, su) with states ^Pi and ^Pi only differ in the sum of the quark 
spins, and can mix via the spin-orbit interaction (and possibly others) to produce 
the physical states i^i(1270) and i^ri(1400). Another example is found in the flavour 
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wavef unctions: the uu, dd and ss are all flavourless (since the flavour of the u cancels 
that of the u) and so there is nothing to prevent the flavour states in a ^5*0 spin/space 
state (for example) from mixing by an annihilation interaction to produce the physical 
states 7r°, rj and rj' . 

By mix, we mean that the basis of quantum number states (e.g. K{}Pi) and 
K(3Pi)) differs from the basis of the physical states (e.g. i^i(1270) and i^i(1400)). 
When we are only talking about two states, we can represent the mixing by a rotation: 



Ki{mQ)- 



( 



\ 



COS Q sin Q 
- sin Q cos d 



(A.5) 



For this particular mixing, we take 9 = 45° for the calculations of Section p.2| . In 
Section we attempt to determine this mixing angle by comparing predictions of 
meson decay models to experimental data. 

In terms of the amplitudes and widths, Eq. |b!].2| becomes with mixing, using a 
decay involving the 7^1(1270)^ for example, 

riVo)+ = ^^I^lf(i270)+I' = cos^ M|^p^)+ +sin^M||3p,)+r- (A.6) 



In Eq. |A.5| , note that we specifled the positively charged mesons, because the 



flavour makes a difference to the mixing angle for these particular mesons. This 
comes about in the quark model Hamiltonian - the term that causes the mixing. 



(if^g ), changes sign when the antiquark instead of the quark is the heavier strange, 
and this leads to a mixing angle of opposite sign when the ^Pi-^Pi Hamiltonian is 
diagonalized (see Eqs. ^]3| and p.4| ). Thus if we have a mixing angle of 45° for the 
fs:i(1270)+ and Ki{UOO)+ (and fs:i(1270)° and Ki{UOO)^), we must use -45° for the 
fs:i(1270)- and 7^1(1400)- (and i?i(1270)° and i?i(1400)0). 
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The other meson mixing that we use in this work^ is due to the tensor part of 
the colour-hyp erfine interaction M: 



I 



is:*(1410) 
fs:*(1680) 



1.00 0.04 
-0.04 1.00 



(A.7) 



^We did initially consider similar mixings of other mesons, but the widths of the decays involving 
these mesons were too small to be included in our results. 



Appendix B 

Meson Wavefunctions Used in this 
Work 



As can be seen from Eq. ^]T|, meson wavefunctions contain four different component 
wavefunctions: space, spin, flavour and colour. The space and flavour wavefunctions 
are given below. The spin wavefunctions are given in Appendix D.3| . The colour 
wavefunction is trivially a colour singlet (see Appendix p.l|) . 



B.l Space Wavefunctions 

In this work we investigate three different sets of meson space wavefunctions, though 
only two are used for calculations. All three are expressed in terms of simple harmonic 
oscillator (SHO) wavefunctions. In momentum-space the SHO wavefunctions are 

i^Tidp) = i^^''{p)yLM,{^,) (B.l) 
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where the radial wavefunctions are given by 

RnL iP) ^ Tin + L + I) [p ) ' ^' ^^-^^ 

L+- 

the oscillator parameter (3 contains all the parameters of the Hamiltonian, and Ln ^ {p^ / 0^) 
is an associated Laguerre polynomial. The quantum numbers n and L are unrelated, 
and have the ranges: n: 0,1,2,..., and L: 0,1,2,... 

In position-space we have the Fourier transform of the above, 

«/,(r1 = i?^r(r)r,M.(f^.), (B.3) 
where the radial wavefunctions are given by 



1 



2 I 



In some cases it may be convenient to pull the p^ or out of the radial wave- 
function, and associate it with the spherical harmonic to create a solid harmonic. 

The first set of meson space wavefunctions that we use for this work are just the 
SHO wavefunction with the same L and Ml as the meson, and the same degree of 
radial excitation (so a meson in its radial ground state, normally denoted with n = 1 
in the n^^^^Lj spectroscopic notation, uses iP^^lml)- '^^^ oscillator parameter [3 is 
taken to be the same for all mesons, and we examine the choice of it in Section ^.1| . 
SHO wavefunctions are expected to be coarse approximations to the true wavefunc- 
tions, since they neglect all but the linear confining part of the meson potential, and 
that they only approximate. However, they are qualitatively similar, and are useful 
for generating analytical results. 

The second set of meson space wavefunctions that we use are those of Godfrey 
and Isgur H, calculated in the highly successful relativized quark model described 
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in Section |1.2.2| . We will label these wavefunctions by RQM. The radial parts of 
these wavefunctions are expressed as linear combinations of the first + 1 radial 
SHO wavefunctions (as in Eq. p.9| ), where was chosen to give a good description 
of the particular wavefunction. The oscillator parameter was fit individually for each 
meson, and the quark masses were also fitted: m„ = 220 MeV, = 220 MeV, and 
m, = 419 MeV. 

The third set of meson space wavefunctions that we use are single SHO wavefunc- 
tions like the first set, but with the effective oscillator parameters (different for each 
meson) of Kokoski and Isgur The effective parameters were obtained by fitting 
the rms momenta of the SHO wavefunctions to be equal to those obtained from the 
second set of wavefunctions. We do not use this third set of wavefunctions for actual 



calculations because of their poor performance in Section R7 



B.2 Flavour Wavefunctions 

We use the following meson flavour wavefunctions in this work: 
For the isovectors: 



7r+ 

For the strange mesons: 

K+ = -us 
K° = -ds 
i?° = -sd 



—ud 
1 

V2 
du. 



uu — dd) 
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K = su. 



For the isoscalars we assumed "ideal mixing" 



0nonstrange \f^^^^ 
'/'strange 

except for the pseudoscalars in a radial ground state, where we assumed "perfect 
mixing" , 

1 1 _ - 1 _ 

^ (0nonstrange Estrange) 1^^^ ~^ ~\J^^ 

f] = ^^^(0nonstrange + Estrange) = — ("UM + (ic?) + — ^SS. 

In terms of the SU(3) flavour representation |l|, 



?78 cos 9p — rji sin 9p 
rjs sin 9p + rji cos 9p 



where 



2 1 

Y 2 *rnonstrange ~l~ estrange 
vnonstrange y g <Pstrange ; 

ideal mixing corresponds to 9p = —54.7°, while perfect mixing corresponds to 9p = 
—9.7°. The latter number is consistent with the angle obtained from the 77-77' mass 
matrix, while data from decay to two photons suggests 9p ~ —20° []T|. The use of 
this second value is considered in Chapter ^. 



Appendix C 

Field Theory Conventions Used 
this Work 

We use the following field theory conventions in this work: 
For the spinor normalizations, 

where if p = {E,p), then p = {E, —p). 
We define the spinor field 

a •> (27r)2 V E 
and use the anticommutators 
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{ba{p),bUp')}^{d^{p),dl{p')} = S'{p-0) 

{ba{p),ba'ip')} = {dc,{p),da'{p')} = 

{bi{p),bUp')} = {di{p),dUp')} = 0, 
which gives for the normahzation of one particle states 

{ba{p)\ba'{p')) = {da{p)\da'{p')) = S'{p-^) 5^o.' ■ 



Appendix D 



Evaluating the Colour, Flavour 
and Spin Overlaps for Models of 
Meson Decay 

In the ^Pq and flux-tube breaking models of meson decay, the overlaps of the colour, 
flavour, spin and space wavefunctions of the mesons and the created pair must be 
calculated. The overlap of the space wavefunctions is accomplished by the integrals 
of Eqs. and p.23| . The other overlaps are discussed below. 

D.l Colour Overlap 

The calculation of the colour overlap is particularly simple, because the mesons and 
the created pair are all colour singlets, with the wavefunction u = -^{RR+GG + BB) 
(using the arbitrary three colours Red, Green, Blue). 

This may be represented in matrix form (this may seem excessive now, but will 
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Ill 



be useful for the flavour overlaps as well) by 



1 ^ 



UJ 



1 

73 



1 
1 



:d.i) 



where the rows indicate a quark of colour R, G, B, and the columns indicate anti- 
quarks of color R, G, B. 

Consider a simple case with only one meson on each side, {uj^2\uj^g). The super- 
scripts indicate that the quark of A (labelled 1) is also the quark of B, and similarly 
for the antiquark (labelled 2). A particular colour combination (e.g. RR) will only 
contribute to the overlap if it is found on both sides of the bra-ket, and the resulting 
term is just the product of the coefficients in the matrices for that colour combination. 
Then the total result is just the sum of these terms over the possible colours, 



E 

Cl,C2 



{^a)ciC2{^b) 



C1C2 — E a) ciC2{^'b) C2C1 ' ' 

Cl,C2 



Tr[cJAt^|], 



:d.2) 



where Ci (02) runs over the quark (antiquark) colours. 

This is easily expanded to the case where there are two mesons on each side (e.g. 



for the first diagram of Figure |2.1|) : 



E (^b)ciC4 (t^c)c3C2 (^a)ciC2 (^^o)c3. 



Cl,C2,C3,C4 

E (^I)c2Ci(^B)ciC4(^ir)c4C3(^c)c3C2 = Tr[^5^BW([wc] 



D.3) 



Cl,C2,C3,C4 



For the case where everything is a colour singlet, the results are 



/ 14 32 1 12 34\ _ / 32 Ml 12 34\ _ - 



32, ,14|, ,12, ,34\ 



(D.4) 
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D.2 Flavour Overlap 



The overlap of the flavour wavefunctions of the mesons and the created pair can be 
calculated using the matrix notation introduced in Appendix p.l| . The flavour of a 
meson may be represented by a matrix where the rows indicate a quark of flavour u, 
d, s,... and the columns indicate antiquarks of color u, d, s,... For example, consider 
the 7r°, 

^10 
0-10 



—=(uu — dd) = —= 
V2 V2 



\ 







(D.5) 



and the K~^, 



—us 



0-1 





(D.6) 



The ^Pq and flux-tube breaking models assume that the pair is created in an 
SU(N) flavour singlet; we take N = 3, since we are not concerned with any mesons 
involving c, b or t quarks: 



1 ^ 



bo = —j={uu + dd + ss) 
v3 



1 

73 



1 
1 



(D.7) 



Taking N to be some other number would change the value of 7 (70) needed to fit the 
data (because of the normalization in Eq. |D.7| ), but would not change the calculated 
decay widths. 

The overlap is given by expressions similar to Eq. p.3| : 



V3 



Trf. 



^A(PB(PC\ 
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Tr[0; 



a9c9o 9b\ 



1 
73 



In many cases only one of these overlaps will be non-zero. 

As an example, consider the flavour overlaps in the decay i^*(892)" 
0i<'*(892)+ = (pK+, giving us 



:d.9) 



K+71 



0. 







(D.IO) 



D.3 Spin Overlap 

Since quarks have spin i, their spins may only combine to give or 1 in a meson. 
They combine in the usual singlet and triplet states (given as e.g. jS'^M^^)): 



|00) 

111) 
|10) 

1 1-1) 



V2 



[[22)12 2) I2 2)122)] 



1 2 2 / 1 2 2 / 
1 

1 2 2 / 1 2 2 / 



1 1 \ 
2 2 / 



) + ll4)IU)] 



(D.ll) 



The overlap of the spin wavefunctions of the mesons and the created pair can be 
calculated by angular momentum algebra. It can be expanded in a sum over possible 
total spins (e.g. for the first diagram of Figure |2.1| ) 



/v^^ Ivl2 34 \ _ 

E {UiM)Sb, U3j2)Sc; SMs\UiJ2)Sa, (jsh)!; SMs) 

S,Ms 

x{SbMs,ScMs^\SMs) {SAMs^l-m\SMs) 



(D.12) 
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where the jVs are all |, but are necessarily written this way to show the connection 
between the left and right sides of the first factor. The first factor is a transformation 
between two different bases for combining four angular momenta, and is independent 
of Ms- On the left side, ji and are added together to obtain Sb, which is added to 
Sc to obtain S" - a similar procedure is carried out on the right side. The 9j symbol 
is defined in terms of this transformation, 



((jlj4)Jl4, (J3J2)J32; JM\{jij2)Jl2, {J^Ja)^^] J M) = 

(_^)2Ji2+2Ji4+i2-i4 + J32-J34^(2Jl2 + l)(2J34 + l)(2Ji4 + 1)(2J32 + 1) 
jl 32 J 12 



(D.13) 



X < 



3a 33 -hi 

JlA J32 J 



which allows us to write the spin overlap in terms of a 9j symbol: 

/v^^ I 12 34 V _ 



-1)^^+73(25^ + 1)(25b + l)(25c + 1) E {SbMsbScMs^ISMs) 

S,Ms 



x{SAMs^l-m\SMs) 



1 

2 


1 
2 


Sa 


1 


1 


1 > 


2 


2 




Sb 


Sc 


S 



(D.14) 



Using an alternative definition of the 9j symbol that couples the quarks differ- 
ently, we can get a similar expression for the spin overlap for the second diagram of 
Figure p]l| which satisfies 



ixfgMsg Xs^cMsc I xfAMs^ Xl\n) 



ix1+5a+Ss+5c /^14 32 ,12 34 v 

^) \XsBMsgXScMs^\XsAMs^Xl-m/- 

(D.15) 
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Eq. p.l5| was used to simplify Eq. 



Appendix E 

Converting to Partial Wave 
Amplitudes for Models of Meson 
Decay 

The decay amplitudes of the ^Pq and flux-tube breaking models derived in Sec- 
tions are given for a particular basis of the final state: 
1^, 0, Mjg, Mj^) = \Q, Mjg, Mj^). Here 9 and are the spherical polar angles of the 
outgoing momentum of meson B (P = Pb) in the CM frame. 

We would prefer to calculate amplitudes for particular outgoing partial waves, 
I J, M, S", L) , since they are what are measured experimentally. Here | J, M) are the 
quantum numbers of the total angular momentum of the final state, \S, Ms) are the 
quantum numbers for the sum of the total angular momenta of B and C, and |L, Ml) 
are the quantum numbers for the orbital angular momentum between B and C. 

The formula for the decay width in terms of partial wave amplitudes is different 
from Eq. |2.5| : 

S,L 
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where 



r^^ = -^|M^^p. (E.2) 



is a partial wave amplitude, and T^^ is the partial width of that partial wave. 
We use two methods to convert our calculated amplitudes to the partial wave 
basis [^]: by use of a recoupling calculation, and by use of the Jacob- Wick Formula. 

E.l Converting by a Recoupling Calculation 

The recoupling transformation directly relates the bases \ fl, Mj^, Mj^) and | J, M, 5", L): 

\J,M,S,L) = Yl (LMlSMsIJM) {JbMj^JcMjJSMs) 

Mjg,M.Jc'Ms,ML 

xjdn YLM.m \n, Mj^,Mj^). (E.3) 
In terms of the amplitudes, the net result we want is the transformation 



M^^(P) = Yl (LMlSMsIJaMj^) (JbMj^JcMjJSMs) 

Mjg,Mjc,Ms,ML 

X JdQ Yl^jJ^Q) M^'a^^^b'^'^Jc (P). (E.4) 

Note that this result holds for any value of Mj^, since M^^{P) is independent of 
Mj^. If we wanted to avoid favouring one particular value, we could sum over Mj^, 
and convert the sum on the left side to a factor of 1/(2 J4 + 1) on the right side. 
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E.2 Converting with the Jacob- Wick Formula 

The Jacob- Wick formula 



I J, M, S,L)=Y1 V ^TTT ^^O-^^l-^^) {JbXbJc -Xc\S\) I J, M, \b, Xc) (E.5) 



relates the partial wave basis | J, M, S, L) to the hehcity basis | J, M, A^, A^), where 
\b and A(7 are the helicities of B and C, respectively. To use it we must first relate 
the helicity basis to the basis Ab, Ac), 



|j,M,AB,Ac) = y^^ JdnDi;;^(<l),e,o)\n,XB,Xc) (e.6) 

_ — * 

where Dj^)^{(j)., 9, 0) is a rotation matrix, and then choose P to he along the positive 
z axis (in the CM frame still), so that A^ = Mjg, Xc = —Mj^, and 1^2, Mjg, Mj^) = 
\fl, Xb, —Xc). 

The final transformation between the amplitudes is 

,/47r(2L + 1) ^ 

^ = 07 {LOS{Mj, + Mj,)\UMj, + Mj,)) (E.7) 

x{JbMj^JcMj^\S{Mj^+Mj^))M^''^^=^^b+^^c)^^b^Jc{Pz). 



Note that Mj^ in the calculated amplitude is replaced by Mj^ + Mj^ . 
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